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Abstract: 

An adequate description of the neutral-current Drell- Yan process at the Tevatron and the 
LHC, in particular, requires the inclusion of electroweak radiative corrections. We extend ear- 
lier work in this direction in various ways. First, we define and numerically compare different 
methods to describe the Z-boson resonance including next-to-leading order electroweak correc- 
tions; moreover, we provide explicit analytical expressions for those. Second, we pay particular 
attention to contributions from 77 and 7-quark collisions, which involve photons in the initial 
state, and work out how their impact can be enhanced by selection cuts. Third, we supplement 
the 0{a) corrections by universal electroweak effects of higher order, such as universal two-loop 
contributions from Aa and Ap, and the leading two-loop corrections in the high-energy Sudakov 
regime as well as multi-photon radiation off muons in the structure-function approach. Finally, 
we present results on the complete next-to-leading order electroweak and QCD corrections within 
the minimal supersymmetric extension of the Standard Model. 
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1 Introduction 



The Drell-Yan-like production of W and Z bosons both provides a standard candle for hadronic 
high-energy coUiders as the Tevatron and the LHC and offers good possibihties to search for extra 
gauge bosons W' and Z' in high-energy tails of distributions (see, e.g., Refs. [ll[2] and references 
therein). For instance, the investigation of the Z-boson resonance, which is well known from 
LEP and SLC experiments, is of great importance for detector calibration, while the analogous 
study of Jacobian peaks of the W boson in appropriate distributions even allow for precision 
measurements of the W-boson mass. Even the effective leptonic weak mixing angle might be 
measurable [3] at the LHC with a precision competing with LEP and SLC. On the theoretical 
side, all these tasks require precise predictions with an inclusion of both strong and electroweak 
radiative corrections and a careful estimate of the remaining theoretical uncertainties. 

The largest corrections are due to strong interactions, mainly described by perturbative 
QCD. The QCD corrections are known to two loops, i.e. next-to-next-to-leading order (NNLO) 
for integrated cross sections [4J and for differential distributions Including corrections up to 
N^LO in the soft-plus-virtual approximation [6] the remaining theoretical error from QCD for 
inclusive cross sections is at the per-cent level or lower. The next-to-leading-order (NLO) QCD 
corrections have been matched with parton showers [7] and combined with a summation of soft 
gluon radiation [8j. 

While QCD corrections to on- or off-shell W- and Z-boson production with leptonic decays 
are very similar, electroweak corrections to the different gauge-boson production processes differ 
considerably. At NLO the electroweak corrections are completely known, both for charged- 
current (CC) [9Vtl2j and neutral-current (NC) [13Vtl7j processes. A tuned comparison of cross 
sections and differential distributions has shown good agreement between the various calcula- 
tions p]ll8t[T9]. Since collinear singularities from photonic initial-state radiation are absorbed 
into the parton distribution functions (PDF), similar to the usual QCD factorization, a pho- 
ton PDF delivers another source of real electroweak corrections. Corrections due to 75 and 
75 collisions arise both in the CC case (W production) [20H22] and in the NC case (dilepton 
production) [MIET]. In the NC case even a leading-order (LO) contribution is induced by 77 
collisions [T6]. Finally, the NLO calculations to the CC Drell-Yan process have been generalized 
to the supersymmetric extension of the Standard Model (MSSM) in Ref. |22j . 

Beyond NLO electroweak corrections, multi-photon final-state radiation has been considered 
both for W-boson [22ll23j and Z-boson production [24]; more recently even multi-photon radiation 
off all charged particles has been matched with the 0{a) corrections in the HORACE program 
in the CC [12j and NC [16j cases. Moreover, the impact of the leading higher-order effects due to 
Aa and Ap as well as the leading two-loop corrections in the high-energy Sudakov regime have 
been investigated for the CC case in Ref. [22] . 

A proper combination of QCD and electroweak corrections is in progress by various groups. 
Different procedures for this combination based on factorization or addition, as implemented in 
HORACE, are described in Ref. p5]. The results discussed there suggest that non-factorizable 
mixed strong-electroweak corrections, which start at the two-loop level, are required in order 
to achieve per-cent accuracy in the predictions. For on-shell Z production part of these 0{aas) 
effects have been calculated in Ref. ^26j . 

In this paper, we complete and extend the existing results on radiative corrections to the NC 
Drell-Yan process in various respects: 

1. We rederive the 0{a) electroweak corrections and document the analytical results for the 
one-loop corrections explicitly. Moreover, we define and numerically compare different 
treatments of the Z-boson resonance in the presence of weak corrections. Specifically, 
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we discuss the "complex-mass scheme" [2Z1I1H], the "pole scheme" [MIEQ], and a scheme 
employing a simple factorization into the LO cross section containing the Z resonance and 
a factor for the weak correction. 

2. We consistently include dilepton production processes involving photons in the initial state, 
which proceed via the partonic processes 77 — >• g7 — )■ + q, and 57 — )■ + q. 
We even take into account the known NLO electroweak corrections [31] to the process 
77 — )• l^l^, which contributes to the LO signal process. 

3. Beyond NLO we consider universal two-loop contributions from Aa and Ap, the leading 
two-loop corrections in the high-energy Sudakov regime, and multi-photon radiation off 
muons in the structure- function approach [32] . 

4. Finally, we calculate the NLO electroweak and QCD corrections within the MSSM. 

For the Standard Model (SM) the presentation in this paper widely follows Refs. |lip22j. where 
the electroweak NLO corrections and the same type of effects beyond NLO are discussed for the 
CC Drell-Yan process. Similarly our discussion of the NLO corrections in the MSSM, presented 
here for the NC case, proceeds along the same lines as in Ref. [22j for the CC case. 

The paper is organized as follows. In Section [2] we set up our conventions and give the 
lowest-order cross sections. Furthermore we describe and discuss the different treatments of the 
Z-boson resonance and the different input-parameter schemes considered in this paper, as far as 
it is necessary for the LO process. In Section [3] the electroweak radiative corrections of points 
1.-3. given above as well as NLO QCD corrections are discussed. The NLO corrections within 
the MSSM are described in Section HI Our discussion of numerical results, which is presented 
in Section [5l comprises integrated cross sections as well as differential distributions for the LHC 
and integrated results for the Tevatron. We also compare our results to results previously given 
in the literature and discuss how effects of incoming photons can be enhanced. Our conclusions 
are drawn in Section [H In the Appendix we describe the factorization of QED-like collinear 
singularities into the photon distribution function, give explicit expressions for the vertex and 
box corrections in the SM, and provide details on the considered SUSY scenarios. 

2 Conventions and lowest-order cross sections 

In this section we set up our conventions for the discussion of the various partonic processes 
contributing to the production of a charged lepton pair. Apart from the Drell-Yan-like process 
qq —7- 7/Z — )• l^l^ and its radiative corrections we consider the photon-induced process 77 — t- l^l^ 
and its radiative corrections. Although the latter does not have a Z resonance, it is an irreducible 
background to qq — )• 7/Z and therefore should be included. The (electroweak) NLO 

corrections to 77 — )• have been calculated in Ref. ^31j) and we only briefly review some of 
the results given there. 

The momenta of the incoming particles will be denoted with pi, i = 1,2, and the ones of the 
outgoing particles with kj, j = 1,2,3. Explicitly we assign the external momenta and helicities 
(cjj, Tj, A) according to 

g(Pl,f^l) + ^(P2,CT2) ^ rikl,Ti) + l+{k2,T2) [+7/g(/c3 , A)] , (2.1) 

7(Pi,f^i) + 7(P2,fT2) ^ r{ki,Ti) + l+{k2,T2) [+7(^:3, A)], (2.2) 

where q generically denotes the light up- and down-type quarks, q = d,n, s, c, b, and I denotes 
the charged leptons I = e, n, t. The possible photons or gluons in the final state deliver part 
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Figure 1: Partonic lowest-order diagrams for pp/pp — s- + X. 

of the real radiation contribution to the NLO corrections. The remaining part of the real NLO 
corrections is induced by the crossed processes q'y/g — >• + q and q^y/g — )■ + q. The 
Mandelstam variables are defined by 

s = (^Jl+P2)^ i={pi-kif, n = (pi-/c2)2, sii = {ki + k2f. (2.3) 

We neglect the fermion masses of the light quarks, m^, and of the leptons, m^, whenever possible, 
i.e. we keep these masses only as regulators in the logarithmic mass singularities originating from 
collinear photon emission or exchange. Obviously, we have s = su for the non-radiative processes 
qq —7- and 77 — t- At LO the Feynman diagrams shown in Fig. [T] contribute to the 

scattering amplitudes. For qq — t- the polarized Born amplitude M]^^ can be written as 

M\9^^- = -i 9l,v9JivXv{s)A^^^f];f^"^A^\ (2.4) 

where e is the electric unit charge, g^fv ^'■^ chiral couplings of the fermions / to the vector 
bosons the functions xv{s) describe the propagation of V , and A'^'^ are "standard matrix 
elements" containing the spin information of the fermions. 

The standard matrix element A"'^ for the quark and lepton chiralities, a = ai = —a2 and 
r = Ti = —T2, is defined as 

•^"'^ = [Vq{P2)l^i^aU^{Pl)] [^^i (^1 ) Tm^t l'K^2)] , (2.5) 

with an obvious notation for the Dirac spinors Vq{p2), etc., and the chirality projectors uj± = 
^(1 ±75). Explicitly the A'^'^ are given by 

A^^ = 2u, A^^ = 2i. (2.6) 

For a fermion / with charge Q f and third component ^. of its weak isospin the left- and 
right-handed couplings to F = 7, Z are given by 

9ff^ = -Qf^ 9ffz = -—Qf, 9ffz = — • (2-7) 

The sine and cosine, Sw and Cw, of the weak mixing angle are fixed by the W- and Z-boson 
masses iWw and Mz as described below in more detail. 
The propagator functions are defined by 

X^{s) = 1, xz{s) = -r^, (2.8) 

where the complex quantities 

^ll = Mi- iMzFz , fi^ = M^- iMwFw (2.9) 
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denote the locations of the poles of the Z- and W-boson propagators (with momentum transfer 
p) in the complex plane. The gauge-boson widths Ty enter the propagator denominators 
only after performing the Dyson summation of all insertions of the (imaginary parts of the) 
gauge-boson self-energies, i.e. using the above xzis) already goes beyond the lowest perturbative 
order. It is well known that this unavoidable mixing of perturbative orders jeopardizes the 
gauge invariance of predictions, in particular in the presence of radiative corrections Before 
describing our solutions to this problem, we recall an important feature of the explicit form of 
the propagator function. While we have chosen a constant imaginary part in the denominator 
of Xzis), the frequently used on-shell (OS) renormalization scheme, as for instance defined in 
Ref. [36] , naturally leads to a running width in the denominator. In the approximation of massless 
decay products of the boson V, the OS version of xvis) is 

^""^'^'^^ ^ s - M2_og + iMv,osry,os X OS x e{s) ' ^^'^"^ 

The two versions of xv{s) are formally equivalent in the resonance region if mass and width of 
V are properly translated 
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1 + 



OS ./i + r2,/M2 



(2.11) 

OS 



Since the W and Z masses and widths are usually quoted in the OS scheme, we shall perform 
this translation before our evaluations. For the masses, the impact of this conversion typically 
is Mz^os — Mz w 34MeV and Mw,os — Mw ~ 27MeV. We perform our evaluation in the 
following schemes for treating the Z-boson resonance, where at this point we describe the various 
procedures only as far as relevant in LO and give the details for the corrections in the next 
section: 

• Complex-mass scheme (CMS): The CMS was introduced in Ref. [27J for LO calculations 
and generalized to NLO in Ref. [28j. In this approach the squared W- and Z-boson masses 
are consistently identified with /i^^ and respectively, i.e. with the location of the poles 
of the propagators in the complex p'^ plane. This leads to complex couplings and, in 
particular, a complex weak mixing angle via = 1 — = /i^/;u|. The scheme fully 
respects all relations that follow from gauge invariance (Ward or Slavnov-Taylor identities, 
gauge-parameter cancellation), because the gauge-boson masses are modified only by an 
analytic continuation. Beyond LO the complex masses are introduced directly at the level 
of the Lagrangian by splitting the real bare masses into complex renormalized masses and 
complex counterterms, so that the usual perturbative calculus with Feynman rules and 
counterterms works without modification. In contrast to gauge invariance, unitarity is 
not respected order by order in perturbation theory. However, spurious terms spoiling 
unitarity are of (N)NLO in an (N)LO calculation without any unnatural amplification, 
because unitarity cancellations, which are ruled by gauge invariance, are respected. More 
details of this scheme can also be found in Ref. [33J- 

In the CMS the LO amplitude (j2.4p is, thus, evaluated with complex couplings gfjy and a 
complex Z-boson mass. 



^More details on this issue, specific examples as well as proposed solutions can, e.g., be found in Refs. 
and references therein. 



4 



• Pole scheme (PS): The PS exploits the fact that both the location ^Uy of the V propagator 
pole and its residue in amplitudes are gauge-independent quantities [29l[39]. The idea [291 
[30] is, thus, to first isolate the residue for the considered resonance and subsequently 
to introduce a finite decay width only in the gauge- independent resonant part. If done 
carefully this procedure respects gauge invariance, but it should be kept in mind that the 
resonant part of an amplitude is not uniquely determined by the propagator structure 
alone, but depends on a specific phase-space parameterization and in most cases also on 
the separation of polarization-dependent parts. A "pole approximation" — in contrast to a 
full PS calculation as performed in this paper — results from a resonant amplitude defined 
in the PS upon neglecting non-resonant parts. 

The LO amplitude (12. 4[) with real couplings defined via the usual on-shell relation = 
1 — = M^/M|, but with the complex Z-boson mass in xz(s), represents the result of 
a particular PS variant. The PS operation here first splits off the polarization-dependent 
structure A'^'^ and subsequently introduces the Z-boson width in the resonant part of the 
form factors via — M|) — t- — M| -|- iMzPz) = Xzis)/s, while the non-resonant 
photon part is not changed. 

• Factorization scheme (FS): Many variants of factorizing resonant structures from ampli- 
tudes have been suggested and used in the literature, but they all share the idea to separate 
a simple resonant factor from a potentially complicated amplitude that does not involve 
resonances anymore. In Ref. for instance, the virtual electroweak correction to Drell- 
Yan-like W production was factorized from the resonant LO amplitude, so that the relative 
correction factor did not involve resonance factors anymoreJl 

For the present case of NC dilepton production we start from the LO amplitude (j2.4p with 
real couplings, as in the PS, and define the relative correction factor for the weak (i.e. 
non-photonic) one- loop correction in the strict limit of vanishing gauge-boson widths. 

We can, thus, compare two different versions of LO cross sections for qq — )• 7/Z — )• l^l^: one 
version delivered by the CMS, another by the PS and FS, which coincide in LO. 

The electromagnetic coupling a = e^/(47r) yields an overall factor to the LO predictions. 
Although the electric charge is always defined (renormalized) in the Thomson limit, the value 
for a can be fixed in different input-parameter schemes. We support the following three different 
schemes: 

• a{0)-scheme: The fine-structure constant a(0) and all particle masses define the complete 
input. In this scheme, the relative corrections to the qq — )• 7/Z — )• cross sections 
sensitively depend on the light-quark masses via a In mq terms that enter the charge renor- 
malization. 

• a{Mz)-scheme: The effective electromagnetic coupling a{Mz) and all particle masses define 
the basic input. Tree-level couplings are derived from a{Mz), and the relative corrections 
receive contributions from the quantity Aa(Mz), which accounts for the running of the 
electromagnetic coupling from scale Q = to Q = Mz (induced by light fermions) and 
cancels the corresponding alnniq terms that appear in the corrections to the qq channels 
in the a(0)-scheme. 



^The relative electroweak correction defined in this way involves the W-boson width Fw only in logarithms 
cx aln(s — + iA/wTw), which result from soft-photon exchange. 
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G ^-scheme: 



The Fermi constant G^j and all particle masses define the basic input. Tree- 

'2 /A#2\ 



level couplings are derived from the effective coupling = '\/2G^M^(l — M^/M|)/7r, 
and the relative corrections receive contributions from the quantity Ar |40j , which describes 
the radiative corrections to muon decay. Since Aa(Mz) is contained in Ar, there is no large 
effect on the qq channels induced by the running of the electromagnetic coupling in the 
G^-scheme either. 

Since light-quark masses are perturbatively ill-defined and can only play the role of phenomeno- 
logical fit parameters, the a(Mz)- and G^-schemes are preferable over the a(0)-scheme for the 
qq annihilation processes. More details on the difference of the three schemes are provided in 
the next section, where we deal with electroweak radiative corrections (see also Ref. 

The differential LO cross section da]^^ /dd is easily obtained by squaring the LO matrix 
element M]^^ , 



dCl ) 



12 647r2s ^ 55 ' 



(2.12) 
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+ 2 QqQi Re [{g'tqzatiz + 9„„z9iiz) + kdtqzdiiz + 9qqz9iiz) Xz(s) 



X2i 



+ 



(bggzl \9iiz\ +\9ggz\ \9uz\ )^ + K\9qqz\ Iduzl +\9qqz\ Wuzl \xz{s)\ 



2\ X2 



'sm2 



The explicit factor 1/12 results from the average over the quark spins and colours, and il. is 
the solid angle of the outgoing in the partonic centre-of-mass frame. In Fig. [2] we show the 
integrated partonic LO cross sections o'^^(s) for the different schemes (CMS and PS/FS) to 
treat the finite Z width, as obtained in the G^-scheme. We also show the relative difference 
rpg/FS = ''■'^^ I PS/FS / <5''"^|cMS — 1 of the results obtained in the different schemes, which turns 
out to be at the 0.01 per-cent level. 

For completeness we state the contribution of 77 — )• 
cross section, 

'd<TLo^^ _ 1 1 

^ 4 647r2s 



dall^/dCl, to the LO differential 



(2.13) 



pol 




For details we refer to Ref. [21]. Here we just mention that we consider 77 — )• l^l^ cross sections 
always in the a(0)-scheme, because the natural scale for the coupling of the external photons is 
Q = 0. In fact, using the a(Mz)- or G^-scheme here would result in large corrections containing 
alnmg terms, which should be avoided. 



3 Radiative corrections to the partonic cross sections in the SM 

In this section we discuss the NLO radiative corrections to the partonic subprocesses contributing 
to the hadronic process pp/pp — )• + X. For the main contribution of qq annihilation, many 
issues discussed here are very similar to the case of e'''e~ — t- 7/Z — )• // as measured in the 
LEP and SLD experiments, for which precision calculations have been performed in the last two 
decades (see, e.g., Refs. [4T1[42| and references therein). 
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3.1 Survey of radiative corrections and calculational details 

The electroweak radiative NLO corrections can be divided into photonic and weak corrections. 
The photonic corrections consist of real and virtual corrections that are induced by the emission 
and exchange of an additional photon. Since only electrically neutral gauge bosons are involved 
at LO, the photonic subset of the complete 0{a) electroweak corrections is separately invariant 
under U(l)eimg gauge transformations. For the qq channel this classification is, e.g., discussed 
in Ref. |35j in more detail, for the 77 channel this separation was introduced in Ref. ^31j. For 
qq annihilation the photonic corrections can be further classified into separately U(l)eimg gauge- 
invariant parts. Specifically, the photonic contributions can be split into initial-state corrections, 
final-state corrections, and interference terms, according to their charge proportionality to Q^, 
Qf, and QqQi, respectively. In this sense the photonic corrections to the 77 channel are final- 
state corrections proportional to Qf. The virtual photonic corrections to the qq channel are 
composed of the one-loop photon exchange diagrams shown in Fig. [3^) and the corresponding 
counterterm contributions; the counterparts for 77 scattering can be found in Ref. [31] • The 
real photonic corrections consist of processes with single-photon emission, qq — )■ 7/Z — t- -I-7 
and 77 — )• + 7, and of the processes q/q'y — )• 7/Z — )• + q/q, which deliver a correction 
to both LO processes qq — )• l^l^ and 77 — )• l^l^. On the hadronic level the photon-induced 
processes are, of course, suppressed due to the smallness of the photon PDF, but on the partonic 
level all processes are of the same order 0{a) compared to the LO processes. Since real photons 
effectively couple with a(0) and since virtual and real photonic corrections are intrinsically linked 
to each other, it is natural to identify the relative coupling of the whole photonic correction with 
a(0), independent of the choice of the input-parameter scheme chosen in LO. This means in qq 
annihilation (and in the crossing-related q/q'y scattering) and in the 77 channel we scale the cross 
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Figure 3: Vertex and box diagrams for the electroweak virtual corrections to qq — )• a) 
photonic corrections, b) weak corrections with hght incoming quarks q = u, d, c, s, b, and c) 
additional diagrams for incoming b-quarks, where G stands for would-be Goldstone boson fields. 



section contributions of the photonic corrections with a(0)a^ and a(0)^, respectively, where a 
depends on the input-parameter scheme as discussed in Section [2j 

The weak 0{a) corrections to the qq channel comprise contributions of the transverse parts of 
the photon, the Z, and the 7Z mixing self-energies {'^'^ , S^^, and S^^), of weak corrections to the 
7/Zgg and 7/Z/~/+ vertices, the ZZ and WW box diagrams, and counterterms. The diagrams 
for the vertex and box corrections are shown in Fig. [8)3) for incoming quarks other than b's. 
For incoming b-quarks, the same diagrams as for incoming d- or s-quarks exist, but in diagrams 
with internal W bosons the b-quark turns into its massive iso-spin partner, the top-quark. For 
this reason, in 't Hooft-Feynman gauge there are additional versions of those diagrams in which 
one or two W bosons are replaced by would-be Goldstone bosons; these diagrams are shown in 
Fig. [3]:). Details and explicit results on the weak corrections to the 77 channel can be found in 
Ref. [31j. In our explicit evaluation we scale the relative weak correction with the coupling a as 
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defined in the respective input-parameter scheme, i.e. the cross section contributions of the weak 
corrections scale like and aa(0)^ in the qq and 77 channels, respectively. 

The NLO QCD corrections to qq — t- are easily obtained from the photonic initial-state 
corrections, i.e. by setting the lepton charge Qi to zero within the photonic corrections, and 
replacing a(0) — )• asinRjCp, with Cf = 4/3 and as(/iR) representing the strong coupling 
constant at renormalization scale /xr. For squared amplitudes with an incoming gluon, q/qg ^ 
7/Z — )■ l~l~^+q/q, we omit diagrams with the external photon coupling to / and replace a{0)Qg — t- 
asinujTp/S with Tp = 1/2 in the respective squared amplitudes with an incoming photon instead 
of a gluon. 

For this work we have rederived the NLO corrections to the qq channel with standard methods. 
More precisely, we performed two independent calculations, with results in mutual agreement. 
The one-loop diagrams and amplitudes are generated with FeynArts versions 1.0 [l3] and 
3.2 [H]. The subsequent algebraic reductions to standard forms are done with inhouse Math- 
EMATICA routines in one version and with the help of FormCalc [35] and FeynCalc [46j in 
the other. In this reduction the appearing tensor integrals are reduced to scalar integrals with 
the Passarino-Veltman algorithm [37]. The scalar integrals are evaluated using the methods 
and results of Ref . |48fl50j , where UV divergences are treated in dimensional regularization and 
the soft and collinear singularities are regularized by small fermion masses and an infinitesimal 
photon or gluon mass m^/g. Since the application of the CMS requires complex gauge-boson 
masses, the results of Refs. |48H50] on the loop integrals had to be generalized accordingly]! 
The amplitude of the virtual correction, Adgq^'"^^ , can be expressed in terms of a "form factor" 

/. virt, (TT j>virt, err , /.virt, err , /.virt, err ,1 t t^- j_ j_ Am- 

fqq = /qg- phot + /gg, weak + /gg, QCD times the LO Duac structurc A , 

II ^virt, err / /.virt, err , j^virt, err , j^virt, err \ A<jr fQ '\\ 

■'^'■qq — Ugg,phot + /gg,weak + /gg.QCoJ ' V"^--*-^ 

For self-energy and vertex corrections this is obviously possible, since all external fermions are 
considered in the massless limit. In 7^ 4 space-time dimensions the calculation of box diagrams 
actually leads to combinations of Dirac chains that are not present at LO. However, since the 
box diagrams are UV finite the four-dimensionality of space-time can be used to reduce all 
Dirac structures to the one of A^^^, as explained in App.[B]in more detail. Finally, we have 
rederived the photonic and QCD corrections for massless external fermions, photons, and gluons 
within dimensional regularization by making use of the results of Ref. |53] for translating the 
IR-divergent scalar integrals from mass into dimensional regularization and by employing the 
statements made in the appendix of Ref. |54j on the structure of rational terms of IR origin. The 
results of mass and dimensional regularization for IR divergences are in perfect agreement. 

Details of our calculation of real photonic (or gluonic) corrections are provided in the following 
section, where we present our results on the photonic and QCD corrections. The contributions 
resulting from the factorization of mass-singular initial-state photonic or gluonic corrections are 
also reviewed there. 

We conclude this overview by summarizing the structure of the hadronic cross section includ- 
ing the full NLO corrections, 



u,d.c,s,b 



^In detail, the scalar two- and three-point functions with complex masses can be explicitly found in Refs. [ST] 
and [48], respectively. The IR-divergent four-point integral of Ref. [49] is also valid for an internal complex mass, so 
that only the regular four-point functions, as e.g. given in Ref. [50] for real mass parameters, had to be analytically 
continued to complex masses. General results on singular and regular four-point integrals with complex masses 
will be published elsewhere [52) . 
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+ 



J daLO(xiPi,X2P2)( 
+ I da5|^cD(^lA,X2P2)+ / dalf}:^cDi^^Pi,X2P2 
+ J d<T5|^hot(^lA,X2P2) + J d&lfp^,,{xiPi,X2P2) 
J2 ft\xi,fil) fjt'\x2,fil)(yj dcr'-^^{xiPi,X2P2) + I daff (xiPi,X2P2) 



g=u,d,c,s,b 



+ 



^)(x2,/x|)(/ da:^0(^^p^^^2P2) + / da'^f{xiPuX2P2) 



(7=u,d,c,s,b 



d(Ti^^(2;i 



Pl,X2P2) (l + ^77,phot + '^77,weak 



+ y dai;'^^l^,,ixiPuX2P2) + J da^"^^*phot(xiPl,X2P2) 
+ (/li ^ /12). 



(3.2) 



Here fa'^^\x, fip) are the NLO PDF for finding the parton a with momentum fraction x in the 
hadron hi with momentum Pj at the factorization scale /Uf • The contributions d'f^L^ 



a6,QCD/phot 



result 



from the PDF redefinitions that describe the absorption of collinear initial-state singularities of 
gluonic or photonic origin (see next section). The factors S"^"^-^ represent the virtual corrections 
to the squared LO matrix elements for the ab initial state, 



(LO 



(3.3) 



3.2 Photonic and QCD corrections 

The issue of a gauge-invariant treatment of the photonic and QCD corrections has been discussed 
in Ref. [35] in detail (including even massive fermions). From the arguments given there and 
the discussion above, it is clear that a consistent way of evaluating the photonic and QCD 
corrections is to use the complex Z-boson mass fiz wherever it appears. Since the weak mixing 
angle is derived from the ratio of the W and Z masses, and Mw does not appear elsewhere in 
these corrections, the quantity Cw can be treated as free parameter in the context of photonic 
and QCD corrections, and Sw as well as the couplings gffz are derived from Cw Specifically, 
we set Cw to jUw//"z in the CMS and to My//Mz in the PS and FS; the numerical difference is, 
however, marginal, as expected. 

The virtual photonic corrections can be decomposed into vertex and box contributions. 



•'gg, phot ~ J qq, phot ' J qq, phot \^ ) 



box, fJT , 



(3.4) 



where the vertex part contains also the photonic contributions to the fermionic wave-function 
corrections. The vertex correction /p^oU'^) consists of an initial- and a final-state part and reads 



/gXphoI(^) = ~y [F'qqV,phot{s) + Fuv,phot{s)\ [QqQl + Q^qZ dlz Xz{s) 



(3.5) 
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Figure 4: Diagrams for real-photon emission. 




q 

7/z 



7 



Figure 5: Diagrams for photon-induced processes with incoming quarks. 



with the renormahzed vertex form factor 



F 



ff V,phot[ 



Qja 
2tt 



+ 2 



(3.6) 

where irrelevant imaginary parts have been discarded. The interference terms of the virtual 
photonic corrections are due to the photonic box diagrams and can be written as 



J qq, photy^' — Jqq \^->^)^Jqq 



{s,i) 



(3.7) 



The correction factors '"^^ are given in App. [Bl 

The real photonic bremsstrahlung corrections to qq — )• 7/Z — )• whose diagrams are 

shown in Figs. |3]and[5l are calculated using the Weyl-van-der-Waerden spinor formalism adopt- 
ing the conventions of Ref. [55]. This results in very compact expressions for the helicity am- 
plitudes M.'^l'^'^'^^'^'^ {X) = J2v='r,z-^ab'v'''^^''^^W^ where ab = qq,qj,qj refers to the partonic 
initial states and V to the exchanged bosons in the respective diagrams. For real-photon emission 
we get 



M 



qqy 



'(+) 



2 V2e%> 9+^(^,2 k2)^ 



Qq 



{ki k2 



Q 



{Pl P2) 



sii- l^v (Pih){p2kz) 

Qq {hk2)* 



sii- f^v {Plh){P2k3) 



s- {ki k3){k2 ks) 

Qi {pi P2)* 

s - /i^ {ki k3){k2 fcs) 



(3.8) 



in the limit of massless fermions, and we have defined /i^ = 0. The spinor products are defined 
by 



i^p • "q 
cos — sm — 
2 2 



e-'^P" cos sin ^ 
2 2 



(pq) = e'^^pAqB = '^VPoQa 
where pA, qA are the associated momentum spinors for the light-like momenta 



Or, 



(3.9) 



p'^ = po(l) sin 9p cos <j)p, sin 9p sin (j)p, cos 6p), 
q^ = qo{l, sin 9 q cos (j)q, sin 6 g sin cos 9 q). 



(3.10) 
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Owing to helicity conservation in the case of massless fermions all amplitudes with ui = 02 or 
T\ = T2 vanish. The remaining six non-zero helicity amplitudes are obtained from the amplitudes 
()3.8p via a parity (P) transformation 



sgn(o-io-2Tir2 



a"l,o-2,ri,r2 



(A) 



and a CP transformation 



9ffV^9ffv 



sgn(aia2nr2)(A^:^;^-^-^^(A)) 



Pl<->P2 



(3.11) 



(3.12) 



Note that in the above formulas the complex masses in the propagators and the couplings are 
not complex conjugated, since P and CP transformations only act on the wave functions and 
momenta entering the amplitudes. 

Apart from the partonic channels with a qq pair in the initial state, we also include the 
photon-induced processes whose diagrams are shown in Fig. [5] for incoming quarks. Of course, 
there are also the corresponding channels for incoming anti-quarks. The amplitudes for the 
photon-induced processes 



9(Pi,cri) +7(P2, A) 
7(pi,A) + q(p2,o-2) 



r{ki,Ti) + l'^{k2,T2) + qik3,a2) 



are related to the ones for real-photon emission by crossing symmetry, 



-A) = -sgn(a2)A^^i''^^'^^'^^(A) 



M 



-5i,o"2,n,''"2 
79 



(-A) = -sgn(o-i)X 



0'l,0-2,Tl,T2 

QQ 



(A) 



P2-0--fc3 



In terms of Weyl-van-der-Waerden spinors the crossing transformation p 
momentum p is obtained by inverting the conjugated parts only, 

PA -PA^ PA^PA- 
The contributions alf'^ and <5"g/|!^ to the partonic cross section are given by 



da; 



-real 



c,ab 



4 2s 



pol 



ab i-^i 



(3.13) 

(3.14) 
(3.15) 
-p of a four- 

(3.16) 
(3.17) 



where the colour factors for the different initial states are Nc^qq = 1/3 and Nc^q-f 
The phase-space integral is defined by 



1. 



d$^ 



d^ki 



d^ko 



(27r)32A;i,o 7 (27r)32A;2,o 7 (27r)32A:3,o 



d^k. 



{27r)^5ipi +p2-ki-k2- ks). (3.18) 



The phase-space integrals in the real corrections (Tg^^' and o'gj'^^ contain logarithmic collinear 
divergences in the limit of massless fermions. Moreover, the real-photon emission integral contains 
a logarithmic soft singularity because of the masslessness of the photon. To regularize the 
soft and collinear singularities we introduce small fermion masses and an infinitesimal photon 
mass according to the generally known factorization properties of the squared amplitudes in 
the singular phase-space regions. This step is usually performed via phase-space slicing, which 
isolates singular regions in phase space, or via a subtraction formalism, which employs an auxiliary 
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function in the whole phase space in order to cancel all singularities. In our calculation we proceed 
as in the treatment of hadronic W production as described in Refs. |lip22). i.e. we employ three 
different methods: (i) soft phase-space slicing with effective collinear factors [11], (ii) two-cutoff 
phase-space slicing |9] for soft and collinear singularities, and (iii) dipole subtraction [561157] . 
Since the detailed formulas for the CC case can be transferred to the present NC case in a 
straightforward way, we do not go into formal details here, but restrict ourselves to the most 
important features of the singularity structure in the final result. 

The analytical results on the photonic corrections to 77 — )• can be found in Ref. [31] . 
The following discussion of final- and initial-state singularities includes both qq and 77 scattering. 

Two types of final-state collinear singularities arise. First, there is a collinear singularity if the 
l^l^ system in the final state receives a small invariant mass Mu, e.g., via a collinear 7* — )• 
splitting. Since we, however, set a lower limit on Mu, this singular configuration is excluded 
from our region of interest. Second, collinear photon radiation off the final-state charged leptons 
is enhanced by the mass-singular factor aln{Q/mi) (with Q denoting a typical hard scale). The 
nature of this singularity is discussed in more detail in Section 13.4.31 where an effective treatment 
of collinear multi-photon emission is described. 

Singularities connected to collinear splittings in the initial state result from q — )• g/79*, q — )• 
s/lQ* (gluon/photon bremsstrahlung in qq annihilation), g/7 — )• qq*, g/7 — )• qq* (gluon/photon 
splittings into qq pairs in g/jq and g/^q scattering), q — )• 57*, q — >■ 57* (forward scattering 
of g or g in jq/q scattering), and 7 — )• l^l^* (photon splitting into pairs in 'yq/q and 77 

scattering). The last splitting corresponds to configurations with a charged lepton lost in 
the beam direction (proton remnant), i.e. it only contributes if not both charged leptons are 
required in the event signature; the contribution of this configuration is enhanced by the factor 
a lii{Q/mi) (again with Q denoting a typical hard scale). The other splittings lead to l^l^ pairs in 
the final state with a gluon, (anti-)quark, or photon lost in the proton remnant; the corresponding 
contributions are enhanced by factors asln{Q/'mq) and aQgln{Q/mq) for gluonic and photonic 
splittings, respectively. These (non-perturbative) singular contributions are absorbed into the 
PDF via factorization, where finite contributions to this PDF redefinition define the factorization 
scheme. In detail the LO PDF fi^\x), describing the emission of parton a out of the hadron h 
with longitudinal momentum fraction x, are split according to (see, e.g., Ref. [9|I58]) 
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ln(4) Pyfi^) - PiM (21n^ + 1) + C^f{z) 



into NLO PDF fi^^ {x, fip), which now include parton emission up to a scale of the order of the 
factorization scale /xp- The splitting functions are given by 



1 - z 



l + {l-zf 



(3.21) 



The coefficient functions Cij{z), defining the finite parts, coincide with the usual definition 
in L)-dimensional regularization for exactly massless partons where the Inniq terms appear as 
1/{D — 4) poles. Details about the photon PDF redefinition are given in App.|Al Following 
standard definitions of QCD, we distinguish the MS and DIS-like schemes which are formally 
defined by the coefficient functions 



r'MS 



Pff{z) ( In 
1 



+ 



9 + 5z 



+ 8z - 1 



DIS/ 



(3.22) 
(3.23) 

(3.24) 
(3.25) 



We use the MRST2004qed PDF [59j which consistently include 0{as) QCD and 0{a) photonic 
corrections. These PDF include a photon distribution function for the proton and thus allow to 
take into account photon-induced partonic processes. As explained in Ref. [58], the consistent 
use of these PDF requires the MS factorization scheme for the QCD, but the DIS scheme for 
the photonic corrections The contributions (^^ab\:iQl)/ phot appearing in (|3.2p result from the 
corrections in the PDF replacements (|3.19p and (j3.20p after these substitutions are made in 
the LO prediction for the hadronic cross section. More precisely, (^^ab'ciCD/phot corresponds to the 
NLO QCD/photonic contribution (i.e. linearized in as or a) proportional to the PDF combination 



X2,)Up) of the partonic ah initial state. 



3.3 Weak corrections 

In the following we sketch the structure of the weak corrections and emphasize those points that 
are relevant for the treatment of the resonance and for the change from one input-parameter 
scheme to another. The correction factor /g^^^^eak' which is introduced in Eq. (13. ID . is decomposed 
according to the splitting into self-energy, vertex, and box diagrams, 

j^virt, fTT j^self, (TT / i vert, err / /vx pbox, err / ^ 2\ (n o/^N 

Jqq,weak ~ Jqq,wea.k\^) + Jqq,weak\^) + /gg.weak V'^; ^ j • [6.Z0) 



*Note that our choice of the factorization scheme and coefRcient functions for incoming photons differs from the 
previously presented results [16i21j . In Ref. [21] the MS scheme was employed, and the coefficient function C^j^(z) 
of Ref. [16) was fixed somewhat ad hoc. Our redefined photon PDF is fixed in such a way that the momentum 
sum rule for the total proton momentum is respected, i.e. our fixation of the photon PDF follows the same logic as 
the gluon PDF redefinition in the DIS scheme for QCD factorization. In the first preprint version of this paper we 
employed a different factorization prescription that also respected the proton momentum sum rule, but was not 
in line with the standard NLO QCD conventions for factorization schemes. Therefore, we switched to the scheme 
described here (see comments at the end of App. 
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The self-energy corrections comprise contributions from the 77, 7Z, and ZZ self-energies, the 
results of which can be found in Ref. [36] in 't Hooft-Feynman gauge. Self-energy corrections 
to the external fermion states are absorbed into vertex counterterms, as usually done in on-shell 
renormalization schemes. The one-loop diagrams for the weak vertex and box corrections are 
shown in Figs. [8)3) and c), their complete expressions are provided in App. [Bl 

(i) Complex-mass scheme 

We first describe the calculation in the CMS. The self-energy corrections explicitly read 



fself.ar (^^ 2 QqQl ^11 ( ^) , ^2 9qqZ dJlZ ^ZZ ( 2 Ql 9qgZ + Qq dflZ ^-/Z f oyN 



where S^^' denote the renormalized (transverse) gauge-boson self-energies, 

^^.■^(s) = j:':p{s) + 6z^^s , 

= J:f{s) - 6 Mi + 6Zzz{s - fil) , 
S^^(s) = S^^(s) + hz^z s + ^SZz^ {s - fil) . (3.28) 

As mentioned above, the explicit results of Ref. j36] on the unrenormalized self-energies 
can be used, however, we stress that complex gauge-boson masses and couplings have to be 
inserted everywhere. The renormalization constants (^M| and 6Zw' are defined in Eqs. (4.9) 
and (4.10) of Ref. [28] for the CMS. They are expressed in terms of gauge-boson self-energies 
and consistently evaluated with complex parameters (but real- valued momenta); in particular, 
no real part is taken in their definition, in contrast to the usual on-shell renormalization scheme, 
as, e.g., defined in Ref. [[36]. 

The vertex corrections can be written as 



(•vert, (TT / ^\ 2 QqQl r™ f"\ , -n^T 2 9qqZ9ilZ \ I -\ , t^r / 

Jqq,we?ik\^l — ' — ^ qq'y ,^eak\S ) + ^ ll-y ,weak\^ I ~^ — 2~ ^ qqZ,^eak\S) + ^ llZ,weak\S) 



(3.29) 



with the renormalized vertex form factors 



-^//y,weak('S) — -^//V,weak(^) + '^//y,weak ' (3.30) 

The explicit expressions for the unrenormalized form factors -F//y,weak('S) are given in App. iBl 
The subscript "weak" indicates that the contributions from photon-exchange diagrams are omit- 
ted both in the form factors and in the vertex counterterms ^'^ffy^eg,\- In the a(0)-scheme, the 
counterterms are given by 



<5//^,weak = - + ^^^11 + -^^/.wcak " 2~Q^^^^^ 

^?/z,weak = ^ + l^Zzz + SZ^^,^^ - ^^SZ.z , (3.31) 



with 



^9ffz = -—Qf [- + ^ — J . ^9ffz = yf- - + — +^9}fz ■ (3.32) 
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Note that the subscript "weak" appears only on the fermionic wave-function renormalization 
constants SZJ , obtained from the the fermion self-energies, because only those receive a photonic 
contribution. We again emphasize the difference between the renormalization constants in the 
CMS |28] and the usual on-shell scheme [36]. In the CMS, all quantities are derived from 
complex masses and couplings, and no real parts are taken from the self-energies that enter 
the renormalization constants. Explicit results can be found in Ref. |28j . In particular, the 
renormalization constant of the weak mixing angle, Ssyf, is connected to the mass renormalization 
of the complex gauge-boson masses. 

The charge renormalization constant 6e/e contains logarithms of the light-fermion masses, 
inducing large corrections proportional to a In m/, which are related to the running of the electro- 
magnetic coupling a{Q) from Q = to a high-energy scale. In order to render these quark-mass 
logarithms meaningful, it is necessary to adjust these masses to the asymptotic tail of the hadronic 
contribution to the vacuum polarization IP'^[Q'^) = Y^T^ {Q"^) / Q"^ of the photon. Using a(Mz), 
as defined in Ref. [60] . as input this adjustment is implicitly incorporated, and the counterterm 
reads 



rCt, cr 



c-ct, a 



a{Mz) 



1 



where 



Aa(Mz) 



TT77 



(0) 



Re{n^^ 



{Ml)} 



a(0) 



a(0) 
~3^ 



Aa(Mz), 



urn 



In 



\my 3 



(3.33) 



(3.34) 



with nj-^j. denoting the photonic vacuum polarization induced by all fermions other than the 
top quark (see also Ref. and Ni = 1 and A'^^ = 3 are the colour factors for leptons and 

quarks, respectively. In contrast to the a(0)-scheme the counterterm S'^ffv\a{Mz) does not involve 
light-quark masses, since all corrections of the form a" ln"'(mj/M|) are absorbed in the LO cross 
section parametrized by a{Mz) = a(0)/[l — Aa(Mz)]. In the G^-scheme, the transition from 
a(0) to Gfj^ is ruled by the quantity AA^\ which is deduced from muon decay. 



V2G„M2 (M| 



a(0)(l + ArW) +0{a^). 



(3.35) 



The counterterm in the G^-scheme reads 



rct, a 
^ffV 



'ffv 



a{0) 



-Ar«, 
2 



(3.36) 



where the one-loop correction AA^^ is evaluated with complex masses and couplings in the 
CMS. This translation of Ar^^^ into the CMS is easily obtained upon analytical continuation 
of the result given in Ref. [36] in the on-shell scheme. Note that Ar^^^ implicitly contains large 
contributions from Aa(Mz) ~ 6% and the (one-loop) correction {c^/s'^)Ap^^^ ~ 3% induced 
by the p-parameter, where Ap^^^ oc G^ra^. Thus, the large fermion-mass logarithms are also 
resummed in the G^-scheme, and the LO cross section in G^-parametrization absorbs large 
universal corrections induced by the p-parameter. In Section 13.4. ll we further elaborate on higher- 
order effects induced by Aa and Ap. 

The box correction fq^^''^^{s, i) is the only virtual correction that depends also on the scatter- 
ing angle, i.e. on the variables t and u = —s — t. The boxes are decomposed into the contributions 
of the ZZ and WW box diagrams. 



box, err / 
,weak v 



/"ZZ, err / ^ 
/go 



+ Jqq [S, t) 



(3.37) 
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The individual correction factors are given in App. |Bl In App. [B] we also give the explicit 
expressions for the vertex and box corrections for incoming b-quarks, where due to the large 
mass of the top quark additional diagrams [see Fig. [St)] have to be taken into account. 

(ii) Pole scheme 

As explained in Section [21 the application of the pole scheme [SHlEn] starts from a fixed- 
order calculation without any special treatment of the resonance. Specifically we calculate the 
weak corrections in the on-shell renormalization scheme of Ref. [36j, i.e. in our PS calculation all 
masses and couplings are real quantities and the Z decay width Fz only appears where it is made 
explicit in the following formulas. The input-parameter schemes are defined as in the previous 
section, with Aa and Ar^^^ derived from real input parameters. 

In a second step, the resonance pole is isolated from the non-resonant remainder and dressed 
by a properly Dyson-summed Breit-Wigner propagator. The definition of a gauge-independent 
residue on resonance, in general, involves some freedom in the more-dimensional phase space, 
because the resonance location fixes only a single invariant. In our case, for instance, two different 
definitions of the residue result if we write the resonant contribution to the LO amplitude either 
as ri{s,t)/{s — M|) or as r2{s,u)/{s — M|) and simply set s — )• M| in the numerators. Off 
resonance (s / M|) the two versions for the residue, ri(M|,i) and r2(M|,u), obviously are not 
the same due to s + i+u = 0. We apply the pole scheme to the form factors /^^'^^eak defined in 
()3.ip . i.e. we single out the resonance pole after splitting off the spin-dependent standard matrix 
elements A'^'^ from the amplitude. Note that the resonant part of /^^"^^eak comprises self-energy 
and vertex corrections only, which merely depend on s, but not on the variables t and u. 

For the vertex corrections this procedure is very simple. The contributions involving Z-boson 

1 /• vert, Z, err 

exchange, Jgg ' ' , 



^vert,Z, err 
•f qq,weak 



is) 



are modified as follows, 

2 9qqz9uz 



qq. 



+ 



Ml 



z9iiz 



-fggZ,weak('5) + PllZ,veak{s) 



^qqZ, wcsiki^z) + -^«Iz,weak(^z) 



^qqZ,wcak(.s) 



- M| + iMzFz 

^<rgZ,wcak(^z) + ^UZ, 



weak 



II Z, weak 



(M 



(3.38) 



while the non-resonant contributions involving photon exchange are kept unchanged. Off reso- 
nance the introduction of the finite Z-decay width Fz in the denominator of the vertex corrections 
changes the amplitude only in 0(0^) relative to LO, i.e. the effect is beyond NLO. 

The treatment of the self-energy corrections is somewhat more involved and requires the 
inclusion of the LO amplitude. The sum of the LO and self-energy contributions is modified as 
follows, 



/< 



LO, o-T 



self, err 
gq,wGak 



Qq Ql 



+ 



9qqZ 9UZ 



M| 



+ 



Qi9lqz + Q>i9liz ^l\s) 
s s- Ml 



QqQl 



V77 



+ 



9qqZ 9UZ 



s-Ml 



\Ml 



Ml 



tf{Ml) 



Ml 
-Ml) 



^^\Ml) 



Ml 
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+ 



Ql 9qqZ + Qq dJlZ 



z) _|_ ('^ 



tf{Ml) 



Qq Ql 



+ 9qqZ 9llZ 



s-Ml 



is) 



tf{Ml) 



{s-Mi)J:^^{Mr 



{s - Mi) 



2^2 



+ 



{Qi 9qqZ + Qq 9llZ 



1 



+ 



S-MI + iMzFz 

,7Z 
'T 



(3.39) 



with E^^(s) = Here we have used the fact that in the on-shell renormahzation 

scheme the renormahzed Z-boson self-energy fulfills ReS^^(M|) = and that the resummed 
terms account for some imaginary parts via Im Sj.^(M|) = MzTz, which holds in 0{a). Off 
resonance the modification changes the amplitude only in O(a^), i.e. beyond NLO. In the 
resonance region (s ~ M^) the terms involving Fz in the denominators do not count as 0{a) 
corrections, but as LO terms. Thus, in order to achieve NLO accuracy there, Fz has to be inserted 
with NLO precision (or better), or the experimental value should be used. Since the residue of 
the propagator is a gauge-independent quantity, this modification can be done in the resonant 
parts without spoiling gauge invariance. In our numerical evaluation we use the experimental 
value. We finally note that the result (j3.39p of the PS substitution can also be obtained upon 
considering the resonance region of an amplitude that results from the full Dyson summation of 
the matrix propagator of the 7/Z system (see, e.g., Refs. [lHHS]). 

The weak box corrections do not become resonant, so that they are not modified in the pole 
scheme. 

(in) Factorization scheme 

As a third option to define the weak corrections, we make use of the fact that the relative 
weak corrections ^I^'^g^j, to the differential partonic cross sections are regular functions of s, 



qq,' 



even in the resonance region (s — )• Mg) without introducing a finite Z width. For the virtual 
photonic corrections this is not the case because of the appearance of corrections proportional 
to aln(s - M|). We, thus, can define the weak NLO correction to the differential partonic cross 
section in the FS scheme by 



da, 



55, weak 



FS 



rvirt 
-'g^jweak 



=0,(5A/z=S|Z(M|) 



(3.40) 



where the LO cross section da^^, as given in (|2.12p . contains the Z resonance structure. The 
subscripts on '^q^'^weak indicate that the Z width is set to zero everywhere and that the Z-mass 
counterterm is derived from the full on-shell Z-boson self-energy (i.e. including both real and 
imaginary parts), in order to avoid double counting of the width effect already present in the 
LO cross section. This simple scheme respects gauge invariance, because the LO contribution 
does (see Section [J]) and the relative correction is derived from the ratio of two gauge-invariant 
quantities, viz. the weak correction and the LO contribution without any Dyson summation. 

As in the PS, the FS calculation only employs real masses and couplings; the width Fz merely 
enters the LO cross section. The input-parameter schemes are defined in complete analogy to 
the PS. 
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Figure 6: Weak corrections ^^^'^^g^k '^77 weak total partonic cross sections for the 

different initial states and the differences Ax between scheme X and the CMS for treating the 
Z resonance. 



In Fig. [6] we show the relative weak corrections (Jl^'^^eak total partonic qq cross sections 



for incoming up- type, down- type, and b-quarks, and (5™Veak ' weak corrections to 77 — )• 

* , is due 

.weak 

quarks. It turns out that the differences 



The difference between 5VI^ , and (5Xr* , is due to diagrams involving W bosons and top 

ddjweak bb,weak o o x- 



Ax = 5:'"' 



gg,weak 



xvirt 

^ ^qq,wcak 



(3.41) 

CMS ^ ^ 



between the different schemes for treating the resonance are below one per mille (apart from 
the W-pair threshold at \/I = 2Mw where only the CMS delivers smooth results). We, there- 
fore, conclude that all three schemes are equally good in describing the weak corrections to the 
Z resonance. 



3.4 Higher-order electroweak effects 

In this section we describe the inclusion of leading higher-order electroweak corrections to the 
parton processes qq — )• 7/Z — )• First, we discuss the inclusion of the leading universal 

higher-order corrections originating from the renormalization of the electroweak couplings in 
the various input-parameter schemes. Second, we consider the leading electroweak corrections 
in the Sudakov regime which are enhanced by large logarithms In^ (s/M-^). Last we discuss 
the inclusion of multi-photon final-state radiation, which is treated using the structure-function 
approach. 
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3.4.1 Leading electroweak effects and choice of couplings 



At moderate scales the leading electroweak non-photonic corrections in the SM are due to the run- 
ning of the electromagnetic coupling = Ana from zero-momentum transfer to the electroweak 
scale, and the large mass splitting between the bottom and the top quark and the associated 
breaking of the weak isospin symmetry. These leading effects are usually quantified by Aa and 
Ap, respectively, and their two-loop effects can be included in a straightforward way, as described 
in Refs. |36 y 58 1 [61 ] l62j. Starting from the calculation within the on-shell renormalization scheme 
with the electromagnetic coupling fixed by a(0), i.e. within the a(0) input parameter scheme 
defined above, the corrections associated with the running of a are included by the resummation 
of Aa via the substitution 

a(0) 



a(0) a(Mz 



(3.42) 



1 - Aa(Mz) 

in the LO prediction, where Aa is defined in (I3.34p . For Ap the leading effects are taken into 
account via the replacements 



+ Ap 



(1 - Ap) c 



(3.43) 



This recipe is correct up to 0{Ap^) and also reproduces correctly terms of C'(Aa(Mz)Ap) [6T1I62] 
in processes with four light external fermions. Note that in 0{Ap^) both one- and two-loop 
corrections to Ap become relevant; explicitly we use the result 



Ap = 3xt [l + (Ml/mf) xt 



1 



2as 



in' + 3) 



3xt 



3V2G, 



mt 



167r2 



Ap 



(1) 



(3.44) 



with the function p^^) given in Eq. (12) of Ref. [63]. In the following we isolate the genuine 
two-loop effects induced by Aa and Ap after properly subtracting the corresponding one-loop 
contributions Aa and Ap(^) already contained in the full NLO electroweak corrections. 

The leading one- and two-loop effects of Aa and Ap in the a(0) scheme are included in the 
LO cross section (j2.12p upon performing the substitutions (j3.42p and (j3.43p . In this context, the 
basic ingredients in (|2.12p are the products a{0)''QgQf, oi{'^)'^QqQiglqz9iiz^ ^.nd oii^Sf [g'^^zdiiz)^ 
of the electroweak couplings defined in (|2.7p . In the following we define g'fjz fo result from g'fjz 
upon applying (j3.43p . Carrying out the above substitutions and keeping terms up to two-loop 
order, the results for these products can be written as 



a(0)2Q2Q2 

ai^fQqQigqqzgJiz 



a{Mzf QlQf = a(0)2 QlQf \l + 2Aa(Mz) + 3Aa(Mz)2 + 

c^iMzf QqQigqnzgJiz 



«(o)' {glqzgJiz? 



ai^YQqQi [gggzaliz [1 + 2Aa(Mz) + 3Aa(Mz)^ 

+ Apa^r + Ap2 + 2Aa(Mz)Apa^r + . . .|, 

«(Mz)=^ irggzgliz? 



«(0)'^ {gggzgJiz? [l + 2Aa(Mz) + 3Aa(Mz)2j + 2Apa''J g'^^^zgJiz 
+ Ap" ((a^r)2 + 26^r gqgzgJiz) + 4Aa(Mz)Ap a^f g'^^^zgJiz 



+ ••• , 



(3.45) 
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where we have introduced the shorthands 



''ql 



r.2 t3 t3 
-w-'w,ij''-'w,Z^ 



+ 



(3.46) 



with YfcT = 2{Qf — fcr) denoting the weak hypercharge of fermion f"^ with chirahty o" = it. 
Dropping the LO contribution and subtracting the relevant one-loop terms, which are propor- 
tional to Aa(Mz) and Ap^^\ the leading two-loop contributions to the coupling combinations 
read 



QqQl 9qgZ9llZ 



{fqqZ9uzf 



a(0) 
hi? 
a(0) 
hi? 



.(1) 



"(0) 



a(0) 
LL2 



a(0)2Q,Q; js^gVftlz Aa(Mz)2 + (^A/, - Ap( 
+ A/52 6-f + 2Aa(Afz)Apa^r|, 

«(0)'{3(5,V5k)' Aa(Mz)2 + 2 (^Ap - Ap^^j"^"^) a^f g^V^k 

+ Ap2 ((a-)' + 26^r5,V5//z) + 4Aa(Mz)Apa^r5,V5k}, 

(3.47) 

where we have indicated the 0(0) input-parameter scheme by superscripts. In Ap*^^^ the super- 
script means which value of a is used in its evaluation. 

The transition from the a(0) to the a(Mz) scheme is rather easy. Since a(Mz) is already 
used as LO coupling, only the replacement ()3.43p applies, but not (j3.42p . Thus, starting from 
the formulas in the a(0) scheme given above, the terms involving Aa(Mz) should be dropped, 
and we obtain for the leading two-loop terms 



QqQl 9gqZ9llZ 
{9qqZ9uzf 



a{Mz) 
a{Mz) 

a(Mz) 
LL2 



0, 



a(Mz)2Q,Qj Ap- ApW 



Q(Afz) 



a^r+ApH^r 



a{Mzfl^2 (^Ap- ApW 



q(A/z) 



'^ql 9qqZ9llZ 



+Ap2 («r)2+26^r9,V5k)}, 



(3.48) 



In the G^-scheme, ac^ effectively involves a factor a(Mz)s^, so that the basic replacements 
read ac^ ^G^s^/ s%f and ffjj^ 9ffZ- This procedure leads to the following leading two-loop 
terms, 



QlQf 



QqQl 9qqZ9llZ 



LL2 



Q^gf{2(^Ap-Ap« 



^ + Ap2^ 



Q^Qz Ap- Ap 



.(1) 



err , ^^W cr r 
'^9' „2 SqqZ9llZ 



I A 2 / 7CTT , ^''W (TT I ''W (T r \ 1 

+ Ap 69/ + agi + (7ggZ5/ZZ ^ , 
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Fi gure 7: Weak corrections (^^g^^gg^j, and "^^^^gg^j^- to the partonic cross sections for the different 
input-parameter schemes, with (dashed hnes) and without leading higher-order corrections due 
to Aq and Ap. 



(3.49) 



We recah that in the CC case [22] the G^-scheme absorbs the full Aa and Ap terms into the LO 
prediction (at least up to two loops), because the CC coupling factor ag^/s^ does not receive 
such universal corrections. In the present NC case this absorption is not complete, and only a 
numerical analysis can assess the size of the remaining explicit universal two- loop corrections. 

In Fig. [7] we show the weak corrections ^uQ^^gak "^df weak partonic cross sections 

for the different input-parameter schemes, including the corresponding higher-order corrections 
due to Aa and Ap. It is clearly visible that the impact of the universal two- loop corrections is 
largest in the a(0)-scheme and smallest in the G^-scheme, as expected. We, therefore, conclude 
that the G^-scheme should be the most stable w.r.t. higher-order electroweak effects among the 
discussed input-parameter schemes. From the above formulas it is also clear that none of the 
schemes is fully optimized to absorb the effects of Aa and Ap into the LO prediction as much 
as possible. While the a(Mz) scheme is more suited for photon exchange, where no leading Ap 
corrections arise, the G^-scheme describes Z-exchange diagrams better, because the generic NC 
coupling e/(swCw) is closer to the weak gauge coupling e/sw than to e. In view of the accuracy 
required for hadron collider physics we do not see, however, the necessity to switch to a mixed 
optimized scheme and take the G^-scheme as default in the following. 

3.4.2 Leading weak corrections in the Sudakov regime 

For dilepton production at large lepton transverse momenta, the parton kinematics is restricted 
to the Sudakov regime, characterized by large Mandelstam parameters s, \t\, \u\ ^ M-^. The 
structure of electroweak corrections beyond 0{a) in this high-energy regime has been investigated 
in detail by several groups in recent years (see e.g. Refs. [641171] and references therein). 

As described for example in Refs. [691171]. the leading electroweak logarithmic corrections, 
which are enhanced by large factors L = ln(s/M^), can be divided into an SU(2)xU(l)- 
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symmetric part, an electromagnetic part, and a subleading part induced by the mass differ- 
ence between W and Z bosons. The last part does not contribute to corrections oc (aL^)" and 
is neglected in the following. The leading (Sudakov) logarithms oc {aL?)"^ of electromagnetic 
origin cancel between virtual and real (soft) bremsstrahlung corrections; for the subleading log- 
arithms such cancellations should strongly depend on the observable under consideration. The 
only source of leading logarithms is, thus, the symmetric electroweak (sew) part, which can be 
characterized by comprising W bosons, Z bosons, and photons of a common mass Mw. In the 
following we consider this type of corrections to the qq annihilation channels of the light quarks, 
i.e. (7 = u, d, c, s, which deliver the dominating contribution to the dilepton cross section. 

The one-loop correction (5g^|sew to the squared amplitude, with chiralities a and r as defined 
above, can be obtained by expanding the full result for the virtual correction 6^^^''^'^ (given in 
the appendix) for large s, \t\, \u\ » M-^. The explicit result can be written as 

Oqq,sew — ^ '-'1,NC -'^'-"1,NC /'-^0,NC/ {<i.OU) 

with factors 

Cq.NC = [gqgZallZ + QqQl) > 

cT^'c^ = {fqqzf + Ql + iahf + Qf + HJ^, 

<NC =|(Co,Nc)'ln(S)+^'^.--5._ln(^) with r 

which have been introduced in Section 8.4.1 of Ref. [71j. In Eq. (I3.50p we did not only include 
the leading Sudakov logarithms oc aL^, but also the related "angular-dependent" contributions 
oc aLln(—t/s) or aLlii(—u/s). Our explicit 0{a) result is in agreement with the general results 
presented in Refs. |69y71j. where the corresponding corrections are also given at the two-loop 
level. These C(a^) corrections can be obtained from the 0{a) result by an appropriate exponen- 
tiation [67j. For the leading "sew" corrections (including a^L^, a^L^ ln(— t/s), and a'^L^ ln{—u/s) 
terms) this exponentiation simply reads [71] 

~ |A^^?|2 exp {4!,ew} = lA^^.^l' (l + + ^lU + • • •) (3-52) 

with ^ 

^gisew = (^-^^ I 2-^^(^rNc)^ ~ -^^C'i^NcC'l,Nc/C'o,Nc|i (3.53) 

where we have suppressed the chirality indices a, r in the notation. 

Particularly in the case of NC fermion-antifermion scattering processes it was observed [70j 
that large cancellations take place between leading and subleading logarithms. In view of this un- 
certainty, we do not include the two-loop high-energy logarithms in our full predictions. Instead, 

(2) 

we evaluate the leading two-loop part 5^^ ,,^^ as a measure of missing electroweak corrections 
beyond 0{a) in the high-energy Sudakov regime. 

Moreover, since the electroweak high-energy logarithmic corrections are associated with vir- 
tual soft and/or collinear weak-boson or photon exchange, they all have counterparts in real 
weak-boson or photon emission processes which can partially cancel (but not completely, see 
Ref. [65]) the large negative corrections. To which extent the cancellation occurs depends on 
the experimental possibilities to separate final states with or without weak bosons or photons. 



tfor4 ,^4_, >0, 
u for I^ji, < 0, 



(3.51) 
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This issue is discussed for example in Refs. [TUITS]. The numerical analysis presented in Ref. [73] 
demonstrates the effect of real weak-boson emission in the distributions in the transverse lepton 
momentum pT,« and in the invariant mass Mu of the lepton pair. At the LHC, at Mu = 2TeV 
the electroweak corrections are reduced from about —11% to —8% by weak-boson emission. At 
PT,i = ITeV the corresponding reduction from about —10% to —3% is somewhat larger. This 
illustrates the sensitivity of weak-boson emission effects to the details of experimental event 
selection, in particular, how dilepton production is separated from diboson production. 



3.4.3 Multi-photon final-state radiation 

The emission of photons collinear to the outgoing charged lepton leads to corrections that are 
enhanced by large logarithms of the form a\ii{mf / Q"^) with Q denoting a characteristic scale 
of the process. The Kinoshita-Lee-Nauenberg (KLN) theorem [71] guarantees that these log- 
arithms cancel if photons collinear to the lepton are treated fully inclusively. However, since 
we apply a phase-space cut on the momentum of the outgoing lepton, contributions enhanced 
by these logarithms survive if the momentum of the bare lepton is considered, i.e. if no photon 
recombination is performed. While the concept of a bare lepton is not realistic for electrons, it 
is phenomenologically relevant for muon final states. 

The first-order logarithm aln{mf /Q"^) is, of course, contained in the full 0{a) correction, so 
that Q is unambiguously fixed in this order. However, it is desirable to control the logarithmically 
enhanced corrections beyond 0{a). This can be done in the so-called structure- function approach 
[32j, where these logarithms are derived from the universal factorization of the related mass 
singularity. The incorporation of the mass-singular logarithms takes the form of a convolution 
integral over the LO cross section a^^, 

aLLFSR= f da^^{pi,p2;ki,k2) dziTlj^{zi,Q^) Qcntizih) f' dz2Tlj^{z2,Q^) ecnt{z2k2), 

J JO JO 

(3.54) 

where the step function Bcut is equal to 1 if the event passes the cut on the rescaled lepton 
momentum Zj/cj and otherwise. The variables Zi are the momentum fractions describing the 
respective lepton energy loss by collinear photon emission. Note that in contrast to the parton- 
shower approaches to photon radiation (see e.g. Refs. [23il24j), the structure-function approach 
neglects the photon momenta transverse to the lepton momentum. 

For the structure function rjf^{z,Q'^) we take into account terms up to 0{a^) improved by 
the well-known exponentiation of the soft-photonic parts [32]; our precise formula can also be 
found in Eq. (2.21) of Ref. |22j . Technically, we add the cross section (j3.54p to the one-loop 
result and subtract the LO and one-loop contributions 



da^^{pi,P2;ki,k2) [ dzi [ dz2 
Jo JO 



6il-zi)6il-Z2) 



+ rj;^'\zi, q2) s{i - Z2) + 5{i - zi) t\I''\z2, q^) 



Qcnt{ziki)Qcnt{z2k2) (3.55) 



contained in ()3.54p in order to avoid double counting. The one-loop contribution to the structure 
function reads 



C'^(^,Q^) = f (\^) (3.56) 



with the variable 

2a(0) 



TT 



m - 1 



In, — ^ 



(3.57) 
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which quantifies the large logarithm. In this context it should be noted that both the full photonic 
one-loop corrections (see Section [3.2p and the multi-photon effects discussed in this section are 
evaluated with a(0) as the photonic coupling in the corrections. Thus, when subtracting the 
one-loop part of Eq. (|3.55p from the full one-loop result, the logarithmic terms oc a{0)\nmi 
cancel exactly in all our considered input-parameter schemes. 

The uncertainty that is connected with the choice of enters in 0{a^), since all 0{a) 
corrections, including constant terms, are taken into account. As default we choose the value 

Q = (3.58) 

with C = 1- III order to quantify the scale uncertainty, we vary ^ between 1/3 and 3. 

4 Radiative corrections to the partonic cross sections in the 
MSSM 

In this section we examine the effect of corrections to pp/pp — )• 7/Z — )• /^/+ -|- X within the 
supersymmetric extension of the SM. A similar study for the case of W-boson production was 
performed in Ref. [2^. Even though Drell-Yan processes do not represent discovery channels for 
supersymmetry it is important to study the influence of SUSY on Drell-Yan processes since they 
will be used at the LHC to calibrate detectors, to monitor luminosity, and to extract information 
on PDF. Measurements on Drell-Yan processes will also allow for precision tests of the SM and 
its extensions through radiative corrections. If there were large corrections due to SUSY particles 
all this information would be biased and therefore not very useful to extract information about 
the underlying physics. 

As an estimate of the impact of supersymmetric extensions of the SM we calculate the SUSY 
corrections to pp/pp — )■ 7/Z — + X within the MSSM. As in Ref. [22] we calculate the full 
MSSM corrections and subtract the SM corrections, so that the MSSM corrections can be added 
to the SM corrections without double counting, 

^gg,SUSY = Sqq,MSSM — 5gg,SM (Afn = M^o) . (4.1) 

Note that we identify the mass Mh of the SM Higgs boson with the mass Mj^o of the lightest 
Higgs boson h'' of the MSSM for the subtraction of the SM corrections. 

We divide the SUSY corrections into the SUSY-QCD and the SUSY electroweak (SUSY- 
EW) corrections. The SUSY-QCD corrections are due to corrections to the qq^/Z vertices as 
shown in Fig. [8] and the quark wave- function renormalization via squark-gluino loops. To obtain 
the SUSY-EW corrections, we proceed as in Ref. [22] and calculate the complete electroweak 
0{a) corrections in the MSSM and subtract the SM corrections. The SUSY-EW corrections can 
be further divided into Two-Higgs-Doublet-Model (THDM) and pure SUSY corrections. The 
THDM corrections are due to the extension of the Higgs sector to two Higgs doublets. It is this 
part of the corrections where we subtract the SM corrections, since in the decoupling limit where 
the mass M^o of the pseudoscalar Higgs boson A*^ becomes large, both the THDM sector of 
the MSSM and the SM Higgs sector coincide if we identify the light neutral Higgs boson of the 
MSSM with the SM Higgs boson. The pure SUSY corrections consist of sfermion, neutralino, 
and chargino loops (see Figs. [9l and [T0|) . 

For the computation of the SUSY corrections we have again performed two independent 
calculations, one using the FeynArts/FormCalc/LoopTools [lHllS] framework and the 
other one using FeynArts and inhouse Mathematica routines. The calculation is done using 
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Figure 8: Example diagram for SUSY-QCD corrections, which are due to squark (g) and gluino 
(g) exchange. 





Figure 9: Example diagrams for pure SUSY verl 
(l), sneutrino (ui), chargino (x), and neutralino 




ex corrections, which involve squark (q), slepton 
(x") exchange. 



X 



X 



Figure 10: Example diagrams for pure SUSY box corrections, which involve squark (q), slepton 
(l), sneutrino (ui), chargino {x), and neutralino (x^) exchange. 



the on-shell scheme as defined in Ref. [36]. Since the LO process is a pure SM process the 
renormalization of Ref. [36] can be applied without modification. To treat the resonance at the 
Z-boson peak we use the LO cross section evaluated in the FS and multiply with a correction 
factor, 

^qq,SVSY = (^qg.SUSY ^^^^ X dg^lps ) (4-2) 

where the relative SUSY correction (5qg^susY) as defined in (j4.ip . can be evaluated without any 
special treatment of the Z-boson resonance, i.e. with a zero Z-boson decay width. We find that 
for the SPS benchmark scenarios [77] (see App. \Q) both the SUSY-QCD and the SUSY-EW 
corrections stay below 2% for partonic centre-of-mass energies up to 2TeV. As an example we 
show in Fig. [11] the partonic LO cross section and radiative corrections for dd initial-states for 
the different MSSM scenarios. 
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Figure 11: SUSY radiative corrections (MSSM — SM) to the partonic process dd — t- Z / 
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5 Numerical results for the cross section pp/pp—)-7/Z^/ 

In this section we describe our numerical setup and discuss the numerical results for the proton- 
(anti-)proton cross section a of the processes pp/pp — )• 7/Z — )• l^l'^ + X in the SM and the 
MSSM. 

5.1 Input parameters and setup for the SM 

The relevant SM input parameters are 





= 1.16637 X 10-^GeV-^ 










a(0) 


= 1/137.03599911, 


a{Mz) 


= 1/128.93, 




= 0.1189, 


-^w,os 


= 80.403 GeV, 


Mz,os 


= 91.1876 GeV, 


Mh 


= 115 GeV, 


rw,os 


= 2.141 GeV, 


Tz.os 


= 2.4952 GeV, 






nie 


= 0.51099892 MeV, 




= 105.658369 MeV, 


rrij- 


= 1.77699 GeV 




= 66 MeV, 




= 1.2 GeV, 


mt 


= 174.2 GeV, 


rud 


= 66 MeV, 


ms 


= 150 MeV, 


mb 


= 4.6 GeV, 



(5.1) 

which essentially follow Ref. [75]. The masses of the light quarks are adjusted to reproduce the 
hadronic contribution to the photonic vacuum polarization of Ref. [76] . The CKM matrix is set 
to unity. We keep finite light-quark masses in closed fermion loops, their numerical impact is, 
however, extremely small in the q(Mz)- and G^-schemes. The 0(a)-improved MRST2004qed 
set of PDF [59] is used throughout. If not stated otherwise, the QCD and QED factorization 
scales are identified and set to the Z-boson mass Mz. 

5.2 Phase-space cuts and event selection 

For the experimental identification of the processes pp/pp — 7/Z — t- l^l^ + X we impose the 
set of phase-space cuts 

Ma > 50 GeV, pt,z± > 25 GeV, |y,±| < 2.5, (5.2) 

where Mu is the invariant mass of the dilepton system, p-p i± are the transverse momenta and yi± 
the rapidities of the respective charged leptons. The cuts are not collinear safe with respect to 
the lepton momenta, so that observables in general receive corrections that involve large lepton- 
mass logarithms of the form aln(m;/Mz). This is due to the fact that photons within a small 
collinear cone around the momenta of the leptons are not treated inclusively, i.e. the cuts assume 
perfect isolation of photons from the leptons. While this is (more or less) achievable for muon 
final states, it is not realistic for electrons. In order to be closer to the experimental situation 
for electrons, the following photon recombination procedure is applied: 

1. Photons with a rapidity \y^\ > 3, which are close to the beams, are considered part of the 
proton remnant and are not recombined with the (anti-)lepton. 

2. If the photon survived the first step, and if the resolution Ri±j = ^{yi± — y-y)^ + 4''i±^ is 
smaller than 0.1 (with </>;±^ denoting the angles between the (anti-)lepton and the photon 
in the transverse plane), then the photon is recombined with the (anti-)lepton, i.e. the 
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momenta of the photon and of the (anti-)lepton are added and associated with the 
momentum of l^, and the photon is discarded. 

3. Finally, all events are discarded in which the resulting momentum of the (anti-)lepton does 
not pass the cuts given in (j5.2p . 

The same recombination procedure was also used in Ref. [22] for single-W production. 

While the electroweak corrections differ for final-state electrons and muons without photon 
recombination, the corrections become universal in the presence of photon recombination, since 
the lepton-mass logarithms cancel in this case, in accordance with the KLN theorem. Numerical 
results are presented for photon recombination and for bare muons. 



5.3 SM predictions for cross sections at the LHC and the Tevatron 

In Tables [1] and [2] we present the integrated LO cross section together with the electroweak and 
QCD correction factors 6ab for the LHC with a centre-of-mass energy ^/s = 14TeV and for the 
Tevatron with y/s = 1.96 TeV. The subscripts a, b of the correction factors 6ab denote the 
various partonic initial states for pp/pp collisions. The results are given for six different ranges 
of the dilepton invariant mass Mn. 

By definition, our LO cross section cr^*-* includes only contributions from qq initial states, 
i.e. we consistently treat all effects from photons in the initial state as corrections. To show the 
effect of the different treatments of the Z-boson resonance (see Section [3.3p we give results for 
the LO cross section in the FS/PS schemes c'^^Ifs/pSj which differ from the CMS only in the 
sub-permille range (< 0.01%). The 77-induced contribution to the LO cross section is given 
separately by the factor 6^^^. Apart from the suppression by the photon PDF, the partonic 
process 77 — t- l^l^ does not involve a Z-boson resonance and therefore is strongly suppressed for 
low invariant mass Mu. However, at the LHC for higher Mu the 77-induced contribution reaches 
up to 5—6% in our default setup. The 0{a) corrections (5^^,phot and ^-y-y^weak have very small 
effect on the integrated cross section at both the LHC and the Tevatron. In Section [5.51 we shall 
pay particular attention to the question how an enhancement of effects of incoming photons may 
be achieved, a question that is interesting for a possible empirical fit of the photon PDF. 

For the photonic corrections we give results for bare muons ((^ggp^ot) and with the recombi- 
nation procedure described in the previous section ((5g|^phot)' where large logarithms oc aln (^jj^^ 
cancel, so that the resulting corrections are smaller. The effect of higher-order final-state brems- 
strahlung beyond 0{a), as described in Section [3.4.31 is small for the integrated cross section, as 
^mui'ti-7 iisver exceeds the 0.1% level. However, as discussed below they become relevant for the 
invariant-mass distribution around the resonance. The correction S'^J[i[_y is given for the central 
scale choice Q = ^/I with an uncertainty estimate obtained from varying the scale Q between 
Q = 3\/I (upper number) and Q = \/I/3 (lower number). Although the g/g7-induced photonic 
processes can be considered as being part of the 0{a) photonic corrections to the gg- induced 
LO process, we do not include them in Sqg^phot, but give them separately by (^g/g-y^phot- They are 
small for all considered Mu ranges in our default setup, as expected from the suppression by a 
factor a and by the photon PDF. 

Our results on the weak correction are given by (5qg,weak- For low Mu the corrections to the 
integrated cross section are of the order of a per cent. For the LHC at high invariant mass 
the weak corrections are enhanced due to large Sudakov logarithms, eventually getting of the 
same order as the photonic and QCD corrections. The smallness of the higher-order weak effects 
i^h.o.weak and the leading two- loop Sudakov logarithms ^gudakov' described in Sections 13.4.11 
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pp ^ l+l- + X at = 14 TeV 



Mii/GeY 


50-cx) 


100-00 
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500-00 
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(To/pb 
<^0|FS/Ps/Pb 


738.733(6) 
738.773(6) 
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0.000303744(1) 
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^77,o/% 


0.17 


1.15 
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-1.81 
-3.34 


—4.71 

-8.85 


—2.92 
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—3.36 
-7.05 


—4.24 
-9.02 


—5.66 
-12.08 
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n /in+0.18 

0-49Io.i5 


n 1 -7+0.06 


n 90+0.07 

^•'^''-0.06 


r, oq+0.09 


n r^+0.13 


Sqq,weak/% 


-0.71 


-1.02 


-0.14 


-2.38 


-5.87 


-11.12 


<^h.o.weak/% 


0.030 


0.012 


-0.23 


-0.29 


-0.31 


-0.32 


4udakov/% 


-0.00046 


-0.0067 


-0.035 


0.23 


1.14 


3.38 


^q/g7,phot/% 


-0.11 


-0.21 


0.38 


1.53 


1.91 


2.34 


'^77^,phot/^ 
'^77,phot/'^ 


-0.0060 
-0.011 


-0.032 
-0.058 


-0.11 
-0.22 


-0.14 
-0.30 


-0.16 
-0.39 


-0.23 
-0.59 


<^77,weaJc/% 
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-0.64 
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4.0(1) 


13.90(6) 


26.10(3) 


21.29(2) 


8.65(1) 


-11.93(1) 



Table 1: Integrated LO cross section and relative correction factors at the LHC for different 
values of the invariant mass cut Mu. 
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pp l+l- +X at ^/s = 1.96 TeV 
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0.00081 
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Table 2: Integrated LO cross section and relative correction factors at the Tevatron for different 
values of the invariant mass cut Mu . 
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and I3.4.2| points towards the stability of the results concerning higher orders in a, especially 
in the resonance region. Following the attitude of Ref. [22] we consider the size of (^sudakov ^ 
a measure for the missing EW higher-order effects beyond NLO. For the LHC this estimate 
indicates a corresponding uncertainty at the level of 1—3% for invariant masses in the range of 

(2) 

1— 2TeV. For the Tevatron (^g^^g^i^Q^ does not even reach the per-cent level up to invariant masses 
of ~ 600 GeV, which suggests that EW effects beyond NLO do not significantly contribute to 
the theoretical uncertainty for Tevatron measurements. 

The NLO QCD corrections (5qcd are evaluated for a fixed scale Hr = fJ-F = Mz and vary 
strongly depending on the size of the cut on the dilepton invariant mass. The statistical error is 
somewhat larger for the QCD corrections since very large cancellations take place between the 
qq and the qg/qg induced channels. 



5.4 SM predictions for distributions at the LHC 

For brevity we restrict our investigation of distributions for the NC Drell-Yan process to the 
situation at the LHC. Already the results for integrated cross sections indicate that the relative 
corrections in the Z resonance region at the LHC and Tevatron are qualitatively very similar. 

Figures [T2 I I13| . and 1141 show differential distributions and correction factors at the LHC. The 
uppermost plots show absolute distributions, followed below by the relative NLO corrections 
(QCD, photonic, and weak corrections) normalized to the qq LO result. Note that for some 
distributions the correction factors are rescaled. The plots at the bottom show the higher-order 
electroweak and photon-induced corrections. The distributions for cr^^^ and a^^'~'''^'^'^ are our 
best estimates for a fJ-^n^ final state and a recombined final state, respectively, and include all 
corrections shown in Figs. [T2t [T3\ and [T31 

In Figure [12] we present the dilepton invariant-mass distribution do" /dMu in the resonance 
and the high- invariant-mass region at the LHC. The Mu distribution on the left shows the well- 
known large corrections due to the final-state photon bremsstrahlung which significantly change 
the shape of the resonance peak, since events belonging to the Z pole are shifted to smaller values 
of the invariant mass Mu . Using photon recombination these corrections are somewhat smaller 
but still in the range of 40—45% at maximum. In the high-energy tail of the distribution, where 
the leptons are an important background for new physics searches, the photonic corrections rise 
in size up to the order of —10%. At 2TeV they are of the same generic size as the QCD and the 
genuinely weak corrections, which are enhanced due to large Sudakov logarithms at high energies. 
At the lower end of the invariant-mass distribution the QCD corrections reach « —170% which 
demonstrates that the inclusion of QCD corrections to NLO only is insufficient there. In order 
to properly describe this end of the distribution, which is determined by the phase-space cuts, 
most probably QCD resummations are necessary — a task that is beyond the scope of this paper. 
The genuinely weak NLO corrections amount to some per cent in the resonance region and 
tend to negative values for increasing Mu, reaching about —10% at Mu = 2TeV. This effect 
is mainly due to the well-known EW Sudakov logarithms. The multi-photon final-state and 
photon-induced corrections around the resonance region are in the range of some per cent and 
thus comparable in size to the weak NLO corrections. In particular, the higher-order multi- 
photon final-state radiation reduces the effect of bremsstrahlung at resonance. The effect of 
universal weak corrections beyond NLO is very small over the whole Mu range. The photon- 
induced corrections are strongly suppressed at the Z pole, but reach the level of a few per cent 
away from the pole. As we observed for the integrated cross section, in the high-energy range 
the 77- induced processes contribute with ~ 5% to o"^^^ in our default setup, where not only 
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Figure 12: Dilepton invariant-mass distribution and correction factors at the LHC in the reso- 
nance region (left) and the high-invariant-mass region (right). 
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Figure 13: Lepton /-"^listribution and correction factors at the LHC in the resonance region 
(left) and the high-px region (right). 
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Figure 14: Dilepton (left) and lepton (right) rapidity distributions and correction factors at the 
LHC. 
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the 77 LO contribution but also the corresponding photonic and weak corrections are included 
in the plots. 

The lepton px-distribution dcr/dpTZ- is shown in Fig. [131 The distribution has the well- 
known Jacobian peak at /- ~ The EW corrections to the /- distribution are similar 
in shape to the CC case [91-112] . The photonic corrections, which are dominated by final-state 
radiation, distort the shape of the peak and are particularly sensitive to the fact whether photons 
are recombined with the outgoing leptons or not. The weak corrections are qualitatively similar 
to the Mil distribution, i.e. they are at the per-cent level and grow negative for increasing 
transverse momenta. As for the dilepton invariant-mass distribution, close to the lower cut on 
P'T,i- the QCD corrections become negative and very large. In contrast to the Mu distribution, 
where the NLO QCD corrections show a moderate size away from the lower end point, the NLO 
QCD corrections to the j?t distribution are insufficient to describe the spectrum at all — an effect 
that is well-known in the literature. The reason for the dramatic rise of the QCD corrections 
for /- ~ Mz/2 lies in the fact that in LO the spectrum receives contributions from resonant 
Z bosons only for i- < Mz/2, but in NLO resonant Z bosons also feed events into the range 
of larger via the recoil of the Z boson against an emitted hard jet in the real corrections. A 
proper description of the transition between the two regions of Z and Z+jet production requires 
careful QCD resummations [8]. At p^z- ~ the NLO QCD corrections are of the order of 
—75% and together with the negative EW corrections this leads to the dip we observe for cr'^'^*"^, 
and at high energies grow huge, reaching the level of several hundred per cent. The photon- 
induced corrections are again small at the peak related to the Z-boson resonance, but reach the 
level of a few per cent at very low and at high p^ . The weak corrections beyond NLO are 
suppressed in the whole considered pT region. The photonic corrections beyond 0{a) reduce the 
size of the NLO photonic corrections for very low /- ^■t the resonance peak, but do not 
exceed 1% in size. 

Figure [m shows the dilepton and the lepton rapidity distributions da/dyu and da/dyi-, 
respectively. The dilepton rapidity yu is defined by 

^"=2^"(4[^j' ^" = ^1 + ^2. (5.3) 

All NLO QCD and EW corrections to the yu distribution are at the level of few per cent, 
indicating the perturbative stability of this observable. This is confirmed by the extremely 
small size (below 0.2%) of the higher-order EW effects shown in the lower left plot of Fig. IHl 
Qualitatively these statements also hold true for the distribution in the rapidity of the lepton, 
where the QCD corrections reach the 5% level. 

5.5 Enhancing effects of incoming photons by cuts 

In this section we study how the effect of photons in the initial state can be further enhanced 
by choosing different phase-space cuts. If the impact of incoming photons can be significantly 
extracted, dilepton production lends itself as a possible candidate for fitting the photon PDF. In 
this discussion, however, it is essential to include also QCD and photonic corrections, which are 
potentially large and especially sensitive to kinematical cuts. We consider the following three 
"77 scenarios" , which are defined by cuts in addition to our default setup, in order to enhance 
the effect of incoming photons: 

(a) Pt,«± < Mil/ 4. 

This cut is motivated by the consideration that 77 fusion involves t- and n-channel di- 
agrams, while the qq annihilation proceeds via s-channel diagrams only, i.e. 77 fusion 
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prefers a small value of sin^*, with 9* denoting the scattering angle of the lepton in the 
partonic centre-of-mass frame, while qq annihilation favours intermediate angles. In LO 
we have Mn = and Px «± = ^VSsinO* = ^MusinO* , i.e. the above cut translates into 
sin0*<i. 

(b) Pt,i± < 50GeV. 

Following the same considerations this cut translates into the condition sin0* < 
100 GeV /Mil for the LO contribution. For increasing invariant masses Mn only smaller 
and smaller scattering angles are included. 

(c) Pt,i± < 100 GeV. 

This case is similar to (b), but with the more relaxed LO condition sin^* < 200 GeV /M;;. 

In Table [3] we present the integrated LO cross section and the photon-induced processes as 
well as the photonic and QCD corrections for the three scenarios. The corresponding dilepton 
invariant-mass distributions and the same types of corrections are shown in Fig. [T5l where the 
bands around the central lines correspond to a variation of the renormalization and factorization 
scale in the range Mz/2 < /iR = /ip < 2Mz. 

The Z-resonance region is clearly dominated by qq annihilation and therefore not suited to 
access effects from photonic initial states. Scenario (a) excludes this region by the applied cut 
completely, since the lower cut on p-p ;± implies Mu > 100 GeV, so that the first two columns of 
numbers in Table [3] are equal. Phenomenologically these two columns are useless, since "edge-of- 
phase-space effects" render QCD corrections unphysically large there. Scenarios (b) and (c) fully 
include the Z peak in the first column of numbers in Table [3l and a significant difference between 
these two scenarios is visible only in the QCD corrections. This is also trivially visible in Fig. [T5l 
The interesting region for studying photonic initial states is the one of higher invariant masses 
Mil, viz. Mil ^ 150 GeV. In scenarios (a) and (b) already for Mu > 150-200 GeV the impact 
of 77 initial states reaches the order of 10-20%; in case (c) it is still 4-7%. In scenarios (a) and 
(b) it should, however, be realized that the effect of 77 initial states is systematically reduced 
by some per cent by the contribution of q/q'j initial states. The correlation between 77 and 
q/qj initial states is also visible in the fact that the sum of the two shows a somewhat smaller 
sensitivity to the variation of the factorization scale than the two individual contributions, as 
seen in the distributions of Fig. [T5j In view of the smallness of the integrated cross sections, 
which are of the order of 0.01 pb to 1 pb depending on the Mu range and the chosen scenario, 
it is not clear whether a fit of the photon PDF within a reasonable accuracy will be feasible. 
This possibility has to be analyzed in a dedicated study that carefully takes into account all 
experimental and theoretical uncertainties, in particular, from higher-order QCD effects. The 
QCD corrections, which are given in NLO in Tableland Fig. 1151 are of the order of —30% to 
—50% and thus quite large; the fact that they are negative, however, helps to further enhance 
the impact of photonic in initial states. The photonic corrections to qq initial states are about 
—10% and thus go into the same direction as well, but the main uncertainty will certainly come 
from QCD corrections. The photonic corrections, which are included in (^^-y o+phot ' weak 
(not explicitly shown here) corrections to the 77 initial states will not play a role in this context. 
In view of the overall size of the cross sections and the sizes of the 77 contributions and the QCD 
corrections, 77 scenario (a) seems the most promising to access the photon initial states. 

5.6 Comparison to SM results of other groups 

In order to make contact to results previously presented in the literature [HIIIIIIIHIIIT] , we have 
added our results on integrated cross sections and NLO EW corrections to the tuned comparison 
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77 scenario (b) 
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Table 3: LO cross section from qq annihilation together with the relative impact from 77 and 
q/qj initial states, as well as from photonic and QCD corrections, evaluated in the three different 
77 scenarios (a), (b), and (c) as described in the text. 
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Figure 15: Dilepton invariant-mass distribution and correction factors at the LHC for the various 
77 scenarios (a), (b), and (c), as described in the text. Within the bands the scale ^ = = /^R 
is varied in the range Mz/2 < fi < 2Mz- 
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0.504(1) 


ZGRAD2 
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LHC, pp — 7- Z, 7 — 7- fi^fj, 




LO [pb] 


NLO [pb] 


6[%] 


HORACE 


739.33(3) 


762.20(3) 


3.09(1) 


SANC 


739.3355(3) 


762.645(3) 


3.1527(4) 


ZGRAD2 


740(1) 


764(1) 


3.2(2) 


our results 


739.343(1) 


762.21(1) 


3.092(1) 



Table 4: Extension of the tuned comparison shown in Table 2 of Ref. [19] for "bare cuts" . 

shown in Table 2 of Ref. [19j. In detail, for this comparison we conformed our input to the setup 
of "bare cuts" described there. Tabled! which shows our results together with the ones obtained 
with the HORACE [16], SANC pT], and ZGRAD2 [HKll] programs, reveals good agreement 
between the various calculations. The remaining differences, which are at the 0.1% level, are 
phenomenologically irrelevant and should be due to slightly different settings in the programs, 
such as the treatment of small fermion masses. 

In order to demonstrate the agreement of our results on distributions with previously pub- 
lished results, a comparison of the genuine NLO EW and multi-photon corrections for various 
distributions to results obtained with HORACE is shown in Figs. [T6l 1171 and [181 In this com- 
parison the HORACE results and the complete numerical setup and input are taken over from 
Ref. [16] . There is very good agreement for the genuine 0{a) corrections, as it should be, because 
these corrections are defined exactly in the same way in the two calculations. Even the multi- 
photon corrections perfectly agree, although HORACE employs a parton-shower approach for 
their modelling, while we use structure functions for collinear multi-photon radiation. The band 
defining our result in Fig. [16] indicates the effect from varying the QED scale (I3.58P by a factor 
3 up and down. A similar comparison in the case of the CC Drell-Yan process was performed in 
Ref. I22], revealing agreement between the two approaches at a somewhat lower level of accuracy. 

5.7 Numerical results on supersymmetric corrections in the MSSM 

Our discussion closely follows the one presented in Ref. [22] for the CC case of single-W pro- 
duction. We choose the SM input parameters and the setup of the calculation (input-parameter 
scheme, PDF, cuts, etc.) as described in Section [5TT] and study the dependence of the corrections 
on the SUSY breaking parameters by considering all the SPS benchmark scenarios |77j . Both for 
the CC and NC case, the generic suppression of the genuine SUSY corrections is insensitive to a 
specific scenario. We therefore refrain from further restricting the SPS coverage by taking into 
account recent experimental bounds in favour of a broader scope in the SUSY parameter space. 
The SPS points are defined by the low-energy SUSY breaking parameters which determine the 
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Figure 16: Comparison of the multi-photon corrections for the dilepton invariant-mass distribu- 
tion to results obtained with HORACE. 
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Figure 17: Comparison of the NLO EW corrections for the dilepton invariant-mass (left) and the 
Pt^i- (right) distributions to results obtained with HORACE, with A showing the difference of 
the two results on the relative correction 5 in the inset, where the HORACE integration error 
defines the error bars. 




Figure 18: As in Fig. [T71 but for the dilepton (left) and lepton rapidity (right) distributions. 
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spectrum and the couplings. For the ten benchmark scenarios under consideration in Ref. [22j 
and in this work, this input [78] is also tabulated in App. ICl 

Dependent SUSY parameters, such as Higgs, chargino, neutralino, or sfermion masses, are 
calculated from the SPS input using tree-level relations. Since the impact of the fermion masses 
of the first two generations is negligible, these masses are set to zero in the calculation of the 
corresponding sfermion mass matrices. Following this approach the SUSY corrections do not 
depend on the lepton generations in the partonic process qq — >• i.e. the SUSY corrections 

presented below are valid both for outgoing electrons and muons. 

In Table [5] we list our results for the SUSY corrections within the MSSM at the LHC. The 
corresponding LO cross sections can be found in Table [TJ We give results for SUSY-QCD and 
SUSY-EW corrections separately as described in Section [H As expected and similar to the CC 
case [22], the corrections for low invariant dilepton mass ranges are negligible at the level of 0.1% 
or below. Only for very high Mu and only for a few scenarios the corrections reach the level of 
1—2%. Similar to the CC case the maximum is reached for the SPS2 scenario where the gauginos 
are particularly light and the squarks and sleptons are so heavy that their negative contribution 
becomes effective only at higher invariant mass. Table [6] essentially shows the same features for 
the situation at Tevatron (for LO numbers see Table [2|), where the SUSY corrections for the 
highest reachable invariant masses are even smaller compared to the LHC. 

In Fig. [19] we show the invariant-mass distribution da/dMu for the LHC. As already observed 
for the integrated cross section, the distributions receive per-cent corrections only for an invariant 
mass in the TeV range. The maximum correction is again found for SPS2 scenario where the 
SUSY-EW corrections reach the 2% level. The SUSY-QCD corrections reach 1% for all but the 
SPS2, SPSS and SPS9 scenarios but never exceed 2% for all scenarios in the considered Mu 
range. 

6 Conclusions 

Neutral-current dilepton production represents one of the most important processes at hadron 
colliders, such as the Tevatron and the LHC. On the one hand, the process acts as a standard 
candle that is, e.g, indispensable for detector calibration and sensitive to the parton distribution 
functions of the proton; on the other, it delivers background to many new-physics searches, such 
as for new Z' bosons. Predictions for this process, thus, ask for the highest possible precision, 
i.e. both QCD and electroweak corrections have to be included as far as possible. 

In this paper we have recalculated and further discussed the next-to-leading order corrections 
in the Standard Model, where we have compared different methods to include radiative correc- 
tions to the Z-boson resonance in a gauge-invariant way. This discussion goes beyond previous 
work, but our numerical results confirm results existing in the literature. The relevant formulas 
are listed explicitly and can be used by other groups. We consistently include channels with 
incoming photons, which starts already with the leading-order contribution of the 77 — t- l^l^ 
process. We also include the electroweak corrections to this process and channels induced by 
IQ/iQ initial states, i.e. all contributions that are formally of electroweak next-to-leading order. 

Beyond next-to-leading order we consider multi-photon radiation off the final-state leptons in 
the structure-function approach, universal electroweak effects, and two-loop electroweak Sudakov 
logarithms at high energies, so that our predictions are of state-of-the-art precision in view of 
electroweak corrections. On the side of the QCD corrections, we include next-to-leading order 
corrections only, so that further improvements via QCD resummations or interfacing parton 
showers are desirable. 
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pp l+r + X at ^ = 14 TeV 





Ma/GeV 


50-00 


lOO-oo 


200-00 


500-00 


lOOO-oo 


2000-00 


SPSla 




-Ew/% 


0.0094(S) 


-0.0041(1) 


-0.053 


-0.43 


-0.33 


0.73 


SPSla 


^SUSY- 


-qcd/% 


0.0060 


0.012 


0.062 


0.34 


1.19 


0.61 


SPSlb 


SsVSY- 


-Ew/% 


0.0076(1) 


0.0021 


-0.018 


-0.31 


-0.67 


-0.20 


SPSlb 


^SVSY- 


-qcd/% 


0.0025 


0.0049 


0.025 


0.13 


0.48 


1.38 


SPS2 


(^SUSY- 


-Ew/% 


-0.046 


-0.12 


-0.40 


0.30 


1.60 


1.94 


SPS2 


i^SUSY- 


-qcd/% 


0.00093 


O.OOIS 


0.0091 


0.045 


0.15 


0.58 


SPSS 


^SVSY- 


~Ew/% 


0.0046(1) 


-0.00072(7) 


-0.021 


-0.32 


-0.66 


-0.20 


SPSS 


^SVSY- 


-qcd/% 


0.0026 


0.0050 


0.026 


0.14 


0.50 


1.37 


SPS4 


^SVSY- 


-Ew/% 


0.013 


0.0005(1) 


-0.061 


-0.24 


-0.24 


0.27 


SPS4 


i^SUSY- 


-qcd/% 


0.0034 


0.0066 


0.035 


0.18 


0.68 


1.30 


SPSS 


(^SUSY- 


-ew/% 


0.025 


0.013 


0.063 


-0.23 


-0.51 


0.11 


SPS5 


(^SUSY- 


-qcd/% 


0.0042 


0.0083 


0.043 


0.23 


0.S5 


1.05 


SPS6 


(^SUSY- 


-ew/% 


0.010 


0.0004(2) 


-0.042 


-0.42 


-0.44 


0.32 


SPS6 


^SUSY- 


-qcd/% 


0.0042 


0.0082 


0.043 


0.23 


0.S5 


1.05 


SPS7 


^SVSY- 


-ew/% 


0.0094(2) 


0.0002(2) 


-0.049 


-0.42 


-0.032 


0.23 


SPS7 


^SVSY- 


-qcd/% 


0.0023 


0.0046 


0.024 


0.12 


0.45 


1.40 


SPSS 


SsVSY' 


-ew/% 


0.012 


0.0039 


-0.035 


-0.34 


-0.17 


0.22 


SPSS 


^SUSY- 


-qcd/% 


0.0017 


0.0033 


0.017 


0.088 


0.31 


1.26 


SPS9 


^SVSY- 


-ew/% 


-0.015 


-0.026 


-0.11 


0.029 


0.095 


-0.0059(1) 


SPS9 


^SVSY- 


-qcd/% 


0.0012 


0.0023 


0.012 


0.059 


0.20 


0.82 



Table 5: Relative SUSY-QCD and SUSY-EW correction factors at the LHC for different values 
of the invariant-mass cut Mu for the SPS benchmark scenarios. The corresponding LO cross 
sections can be found in Table [TJ 
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pp l+l- +X at ^/s = 1.96 TeV 





Mu/Ge\ 


50-00 


lOO-oo 


150-00 


200-00 


400-00 


600-00 


SPSla 




-Ew/% 


0.0022(2) 


0.0034(1) 


0.052 


0.032 


-0.14 


-0.26 


SPSla 


^SVSY- 


^qcd/% 


0.0059 


0.011 


0.029 


0.047 


0.15 


0.30 


SPSlb 


<^SUSY- 


-Ew/% 


0.0043(1) 


0.0050 


0.027 


0.021 


-0.064 


-0.21 


SPSlb 


(^SUSY- 


qcd/% 


0.0025 


0.0046 


0.012 


0.019 


0.059 


0.11 


SPS2 


^SUSY- 


-Ew/% 


-0.051 


-0.12 


-0.30 


-0.40 


0.048 


0.41 


SPS2 


i^SUSY- 


qcd/% 


0.00092 


0.0017 


0.0043 


0.0070 


0.021 


0.040 


SPS3 




-Ew/% 


0.0017(1) 


0.0027 


0.026 


0.018 


-0.074 


-0.20 


SPS3 


^SVSY- 


qcd/% 


0.0026 


0.004S 


0.012 


0.020 


0.061 


0.12 


SPS4 


<^SUSY- 


-Ew/% 


0.0076(1) 


0.0030 


0.015 


-0.010 


-0.21 


-0.059 


SPS4 


i^SUSY- 


-qcd/% 


0.0034 


0.0063 


0.016 


0.026 


0.081 


0.16 


SPS5 


<JsUSY- 


-ew/% 


0.020 


0.024 


0.13 


0.14 


0.037 


-0.0027(1) 


SPS5 


(^SUSY- 


qcd/% 


0.0042 


0.007S 


0.020 


0.033 


0.10 


0.20 


SPS6 


^SUSY- 


-ew/% 


0.0048(1) 


0.0044 


0.034 


0.019 


-0.14 


-0.27 


SPS6 


i^SUSY- 


qcd/% 


0.0042 


0.007S 


0.020 


0.033 


0.10 


0.20 


SPS7 


<^SUSY- 


-ew/% 


0.0054(1) 


0.0040 


0.025 


-0.0017 


-0.21 


-0.30 


SPS7 


(^SUSY- 


qcd/% 


0.0024 


0.0044 


0.011 


0.018 


0.056 


0.11 


SPSS 


5suSY- 


-ew/% 


0.00S4(1) 


0.0067 


0.021 


0.0059 


-0.15 


-0.19 


SPSS 


^SVSY- 


-qcd/% 


0.0017 


0.0032 


0.00S2 


0.013 


0.040 


0.077 


SPS9 


^SUSY- 


-ew/% 


-0.017 


-0.022 


-0.02S 


-0.067 


-0.014 


0.11 


SPS9 


(^SUSY- 


qcd/% 


0.0012 


0.0022 


0.0056 


0.0090 


0.027 


0.052 



Table 6: Relative SUSY-QCD and SUSY-EW correction factors at the Tevatron for different 
values of the invariant-mass cut Mu for the SPS benchmark scenarios. The corresponding LO 
cross sections can be found in Table [2j 
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Figure 19: Dilepton invariant-mass distribution and relative SUSY-QCD and SUSY-EW correc 
tion factors at the LHC for the SPS benchmark scenarios. 
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Our discussion of numerical results comprises a survey of corrections to integrated and differ- 
ential cross sections, which shows the impact of the various types of corrections in detail. In this 
context we pay particular attention to effects from incoming photons, because it is not yet clear 
from results in the literature whether these effects are phenomenologically important or swamped 
by other uncertainties as, e.g., originating from QCD effects. Our results suggest that effects 
from 77 initial states should be significant in the invariant-mass distribution of the dilepton pair 
and that these effects can be enhanced to the level of 10—20% by appropriate cuts. This study 
is, thus, particularly interesting for a possible future fit of the photon distribution function that 
is part of the DGLAP evolution in the presence of electroweak corrections. 

Finally, we have presented results on radiative corrections within the minimal supersymmetric 
extension of the Standard Model, which have not yet been known in the literature. Similar to 
known results on charged-current W-boson production, the difference to the pure Standard Model 
corrections is small for viable supersymmetry scenarios. In detail, the supersymmetry corrections 
are below the 0.1% level and thus irrelevant near the Z-boson resonance and grow only to about 
1—2% in the TeV range for the dilepton invariant mass. These results confirm the role of dilepton 
production as a standard candle at the Tevatron and the LHC. 
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A Factorization of the photon PDF 

Since the issue of redefining the photon PDF in order to absorb the collinear divergence arising 
from the splitting q — >• 57* of incoming quarks (or anti-quarks) seems not to be completely 
settled in the literature yet, we here give a brief account on the analytical structure of this 
collinear singularity both in dimensional and mass regularization and derive the corresponding 
factorization formula for the photon PDF in the MS and DIS factorization schemes. To this end, 
we can make use of existing results, e.g., available in papers on dipole subtraction to deal with 
this collinear singularity. 

In mass regularization this type of singularity is treated in Section 5 of Ref. [57] in detail. 
After extracting the singularity from the full phase-space integral by simple subtraction, the 
singular integration is performed analytically with a small quark mass ruq and added back to the 
result for the cross section. By construction this procedure can be performed with a "spectator" 
particle in the initial or final state, where its choice is arbitrary for this type of splitting. Writing 
the (unpolarized) partonic cross-section part that is to be re-added generically as 

'^Ta-^qxiPq^Pa) = dx 'H{Q'^ , x) a^a^x {k^ = XPq,Pa), (A.l) 

where a^a^x is the partonic cross section of the related process with an incoming photon instead 
of the q — )• g7* splitting. The momenta pq and Pa correspond to the incoming quark q and the 
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other massless particle o, respectively, and the functions % read (cf. Eqs. (5.19) and (5.27) of 
Ref. [57]) 



^W'«(s,x) = In 



s(l - xf 



1 — X 



nf[P\x) = ln( - 



S = {Pq+Paf, (A.2) 
P2 = (p^. _ /,^)2, (A.3) 



where refers to the situation with the incoming particle a as spectator and T-L^ case 

with a massless final-state spectator j of momentum pj. 

The respective results in dimensional regularization can be deduced from the original Catani- 
Seymour paper [79], where the dipole subtraction method was introduced for massless QCD. The 
case with an initial-state spectator is described in Section 5.5 there, and the case with a final- 
state spectator in Section 5.3. Starting from Eq. (5.152) in the former case and from Eq. (5.73) 
in the latter, and identifying the QCD coupling factor CpOs with the QED factor Q^a, we find 



2\ e r 



r(i-e)V s 

^ '^^ r(l-e)^-P2 



P^f{x) + 2P^f{x) ln(l -x) + x 
P^f (x) + P^f {x) In i- I + X 



(A.4) 
(A.5) 



in = 4 — 2e dimensions. Note that in the last relation we had to translate the kinemati- 
cal variable (pkPa) of Ref. ^79j into our kinematics, which is accomplished by the replacement 
(PkPa) iPjPq) = {Pjk'y)/x = -P^/(2x). Comparing Eq. ()A.2p with Eq. ()A.4p . or Eq. ()A.3P 
with Eq. (1X5]) . we see the following correspondence between the collinear divergence in mass 
and dimensional regularization. 



(4V)^ 1 

fa 



-p 



e e 



7/1 



o (lnm^ + 21nx + l)P^/(x), 



(A.6) 



which is obtained in either case, i.e. with a spectator in the initial or final state, as it should 
be. The correspondence ()A.6P can be used to translate a result for the collinear singularity of 
the initial-state splitting q — )• q^* from dimensional to mass regularization, or vice versa, in any 
procedure to treat the collinear divergence, i.e. it is universally valid also in other subtraction 
procedures or in methods employing phase-space slicing. Applying this correspondence to the 
photon PDF redefinition in dimensional regularization. 




(4vr) 



r(i-e) 



P,f{z)+C^f{z) 



(A.7) 

we obtain the result (|3.20p for this redefinition in the mass regularization scheme, where the co- 
efficient function C^j (z) is given by Eq. (j3.25p for the the two considered factorization schemes 
(FS), independent of the chosen regularization. Since the MS scheme merely rearranges the 
divergent terms in dimensional regularization (plus some trivial universal constants), the coeffi- 
cient function obviously vanishes, C^{z) = 0. The fixation of Cj^j^{z) in the DIS scheme is less 
trivial. Similar to the gluon PDF redefinition in the DIS scheme for NLO QCD corrections, the 
choice p.25p ensures that the whole PDF redefinition does not change the proton momentum 
(see e.g. Eq. (6.9) of Ref. [79]). 
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Finally, we mention a subtle point in the fixation of the coefficient functions in the DIS 
scheme. The DIS factorization scheme is defined in such a way that the DIS structure function 
F2 from electron-proton scattering does not receive any corrections. This condition uniquely 
fixes the (anti-)quark PDF redefinition. The gluon PDF redefinition is performed in such a way 
that the total proton momentum remains unchanged after the PDF redefinition. The simplest 
choice obviously is to subtract the same z-dependent finite part from the gluon PDF that was 
added to the (anti-)quark PDFs; this is expressed by the analogon of Eq. (j3.25p in NLO QCD. 
In NLO QED the role of the gluon is taken over by the photon, which finally leads to Eq. (j3.25p . 
Obviously, employing the sum rule alone as criterion to fix C}^j^{z) is not unique, but involves 
some convention, since an integral over z does not fix the z dependence. A different choice would 
result in mass regularization if the coefficient functions C^^^{z) were just defined to quantify 
the finite parts in addition to the divergent contributions that are proportional to InrrigPabiz). 
The result for the redefined photon PDF would still satisfy the proton momentum sum rule, but 
the z dependence of the cross section with the incoming photon would be somewhat different. 
In the first preprint version of this paper such a choice was in fact made, and we found results 
on the contributions from ^7 scattering that differ from the results of this paper at the level 
of up to a few per cent in the off-shell tail of the Z-boson resonance; for observables that are 
dominated by resonant Z bosons effects due to q'j scattering are negligible in either scheme. The 
"correct" choice, i.e. the one that is in line with the standard definitions made in NLO QCD, 
is obtained upon first translating the mass-regularized divergence from the q — >• q^* splitting 
into dimensional regularization via the correspondence (jA.6|) and defining the remaining part as 
coefficient functions, as expressed in Eq. ()3.20p . 

B Vertex and box corrections 

In this section we state explicitly the expressions for the vertex form factors -^Jy-yweak' ^ = 7) ^ 
(j3.29p . and the photonic (|3.7p and weak box diagrams (j3.37p . The occuring scalar integrals Bq, 
Co, and Dq depend on their arguments as follows, 

Bo{pl,mo,mi) = ^'^'^^\ / d^g o , -mr/ ^, ^ 2 , -m ' (^-1) 



ivr^ J [g2 _ 77^,2 _|_ jQ] Jl'g _|_ _ jy{2j_ io] ' 

Co{pI, {P2 -pif,pl,mo,mi,m2) = ^^—^ [ d^q 

^ [g2 - m2 + iO] [(g + pi)2 - mf + iO] [{q + pa)^ - + iO] ' ^^''^^ 



Do{pI, {P2 -Plf, (P3 -P2f,pl,pl, iP3 - Plf ,7710,7711,1712,1713) = ^^^^^ [ d^q 

1 

^ 

[g2 -ml + iO] [(g + pi)2 - + iO] [{q ^ P2Y - iO] [(g ps)^ - m-l + iO] ' 



(B.3) 



with D denoting the number of space-time dimensions. The evaluation of these scalar integrals 
with real or complex masses has been briefly described in Section [3j For the form factors for the 
weak vertex corrections we obtain 
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(B.4) 



F 



a 



//7,weak 



2g/(2s + M^)So(0,0,^w) 



Sir [sslfQf 

+ (2llj -Qf) (2s + {3s + 2n^)Bo{s, 0, 0) + 2{s + /x^)'Co(0, 0, s, 0, /xw, 0)) 
- 2llj ((s + 2hIj)Bo{s, /xw, /^w) - 2/xiv(2s + /x^)Co(0, 0, s 



2 \2 



+ 



-4(s + /x|)2Co(0,0,s,0,/xz,0) 



4s + 4(2s + /x|)Bo(0, 0, Mz) - 2(3s + 2/x|)5o(s, 0, 0) 



(B.5) 



//Z,weak 



(^) = +^ 



2(4,/ - Q/4)(2^ + /^w)^o(0, 0, Mw) 



+ {Iljc^ - Iljsl + Q/4)(2s + (3s + 2/x^)Bo(l, 0, 0) 
+ 2(s + ^^)2Co(0,0,s,0,^w,0)) 

j((s + 2/i^)Bo(s, /xw, /^w) - 2/x^(2s + Mw)Co(0, 0, s, /^w, 0, /^w)) 



+ 



(4,/ - Qfs 



2 \2 



4s + 4(2s + m|)^o(0, 0, Mz) - 2(3s + 2//|)5o(s, 0, 0) 



4(s + ;u^)^Co(0,0,s,0,Mz,0) 



(B.6) 



Due to the heavy mass of the top quark the weak corrections for incoming b-quarks differ 
from those for partonic processes with incoming d/s- quarks by <^-^bb7/z weak ~ -^bb7/Zweak ~ 
^ddM7/z,weak' ^here 



bb7,weak 



(^) ='^^bWcak('?)=0, 



(B.7) 



'^-^bb7,weak(«) = +^ 



a 



(2S + /xw)(Bo(0,0,Aiw) - -Bo(0,mt,/iw)) 

+ (3s + 2m^)(So(s, 0, 0) - 5o(s, mt, mt)) 

+ 2(1 + //|^)2(Co(0, 0, s, 0, Mw, 0) - Co(0, 0, s, mt, mw, "Zt)) 

+ 3//w(2s + Mw)(Co(0, 0, s, //w, 0, yuw) - (^0(0, 0, s, //w, "it, Mw)) 



(mt +/Xw)(-Bo(0,mt,/xw) + 2-Bo(s,mt,mt) - 3Bo(s, /xw, /^w)) 



3 5 
+ s(So(s,mt,mt) - -So(s,|Uw,A*w) " 2^ (Mw(2s + 3/Xw) - "it("^t +«)) 



bbZ,weaJj 



(^) = +^ 



X (2Co(0,0, s,mt,Atw,"^t) + 3Co(0, 0, s, //Wj^^t, A*w)) 
1 



(B.8) 



2(2s + /x^)(2s4 - 3)(So(0, 0, /xw) - So(0, mt, /xw)) 



87r[ (2s2,-3)ss2, 
+ (3s + 2/x^)(4s^ - 3)(So(s", 0, 0) - Bq{s, mt, mt)) 
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- 12c4/Ltw(2s + /iw)(Co(0,0,s,//w,0,/xw) - Co(0, 0, s, //w, "it, ^w)) 
+ 2(44 - 3)(s + //^)2(Co(0, 0, s, 0, /xw, 0) - Co(0, 0, s, mt, /^w, "^t)) 



2 



(m? + ^w)(4^o(s,?Tit,"it)si - 3Bo(s,m,/^w)('S? 



(2sw 3)ss^^|y 

+ 5o(0,mt,^w)(2si - 3)) + (2ss|, - 6/iw)-Bo(s, mt, mt) + •sQ - 5sw^ 

+ - - Cw)^-Bo(s,/xw,/Uw) - 12s/UwCo(0,0,s,mt,/iw,"it) 

- 3 - '2ml fi^ - "^t •5) (Co(0, 0, s, mt, ^w, "^t) + (^0(0, 0, s, /iw, "T-t, Aiw)) 

- 3 (3/iw - "it) ^"0(0, 0, s, fiw,mt,fiw) + 2s^(/Uw(2s + 3^^) " "^t ("^t + «)) 



X (2Co(0, 0, s, mt, /iw, "it) + 3Co(0, 0, s, /xw, "T-t, Mw)) 



(B.9) 



As mentioned in Section [3.H the calculation of the box diagrams leads to additional Dirac 
chains. However, these can be reduced to the Dirac structure (|2.5p appearing in the LO matrix 
element, 

A"^ = [vgj'^u^ug] [uaf^co^vi] . (B.IO) 

We used the identities [2H] 

[VqftiUJ±Uq][ui^2i^±Vl] = +uA^^/2, 

[vq^ioj±Uq][ui^2u;^vi] = (B.ll) 

which are valid in four space-time dimensions, to reduce the Dirac structures of the box diagrams. 
The electroweak corrections due to box diagrams can be written in terms of the functions 

bf'^{Mv,Mv',mQ,mq,mi) = -2^Co{mf ,mf , s, Mv,mi, My) + Co(m^, m^, s. My , mg. My/) 

- {i - rnQ)Do{ml,mg,mf ,mf , s,i, Mv,mQ, Mv' ,mi)^ , (B.12) 

bt^{Mv,Mv',mQ,mq,mi) = ^{2u{Bo{s, My, Mv) - Bo{i,mQ,0)) 

- t{ml - M^ - M^, - i + u){Coimf,ml,i, 0, My.mQ) + Co(mf , m^, £, 0, My,, wq)) 

- (t^ + + s(m^ -M^ -M^,))(Co(mf,mf,s,My,mi,My/) 
+ CQ{ml,ml, s, Mv,mQ, My/)) + (t{M^ + M^, - - 2sf 

+ {u{2u - m^) - 2s^){M^ + M^, - m^Q - 2s) - 2u{iP' - M^M^,) + us{s - m\) - s^) 
X Do{m'^,rn'^,m'f ,7711 , s,i, My, ruQ, Mv' ,mi)^ , (B.13) 

bt^{Mv,Mv',mQ,mq,mi) = -6^^(My , My/, mg, m^, mOI^^^ , (B.14) 
bt^{Mv,Mv',mQ,mq,mi) = -bf^{Mv,Mv',mQ,mq,mi)\f^^. (B.15) 

Here the subscripts q and / refer to light quarks and leptons. The heavy mass mg of the weak 
isospin partner of down-type quarks is present only in the case of WW box diagrams for incoming 
b-quarks. In terms of these functions the electroweak box correction factors read 

fll^^is^i) = a2(5^-^^5[,^)2(5-(^z,/xz,0,0,0) + 6S^(/iz,/xz, 0,0,0)) (B.16) 
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for the ZZ box-diagrams, and 



J. WW, +- 
Jqq 






) = o, 


rWW, -- 
Juu 


-{sj) 


= 4,4 ^« (MW 


/iw, 0,0,0) , 


WW, — 
•'dd 


-{s,i) 




//w, 0,0,0) , 


^WW, -- 
•'bb 




= 4 4 &i (/^W 


//w, "It, 0,0) 



(B.17) 
(B.18) 

(B.19) 

(B.20) 



for the WW box-diagrams. The photonic box corrections are given by 



fqq'"'' = o^QgQf {bt'^{m^,m^,0,mq,rni) + bl''{m^,m^,0,mq,mi)) , (B.21) 

fqq'"^ = '^O? QqQlfqqZduZ {b^ {fJ'Z, m^, 0, rUq, TTll) + b^'^ {flz, m^, 0, rUq, TTll)) . (B.22) 

C SPS benchmark scenarios 

For the SPS benchmark [77] scenarios discussed in this work we use the low-energy input specified 
in Table [71 The input variables are the ratio tp of the vacuum expectation values of the Higgs 
bosons giving rise to up- and down-type fermion masses, the mass of the CP-odd Higgs boson, 
M^o, the super symmetric Higgs mass parameter /i, the electroweak gaugino mass parameters 
Ml 2, the gluino mass rrig, the trilinear couplings Ar,t,b; the scale at which the DR- input values 
are defined, /iR(DR), the soft SUSY-breaking parameters in the diagonal entries of the squark 
and slepton mass matrices of the first and second generations Mfi (where i = L,R refers to the 
left- and right-handed sfermions, f = q,l to quarks and leptons, and / = u, d, e to up and down 
quarks and electrons, respectively), and the analogous soft SUSY-breaking parameters for the 
third generation M|p. 
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